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Abstract. The main purpose of this paper is to give a survey of recent developments on a 
^ ^ classification of soliton solutions of the KP equation. The paper is self-contained, and we give 

5-H complete proofs of theorems needed for the classification. The classification is based on the totally 

non-negative cells in the Schubert decomposition of the real Grassmann manifold, Gr(A'^, M), the 
set of Af-dimensional subspaces in M^^. Each soliton solution defined on Gt{N, M) asymptotically 
consists of the N number of line-solitons for 3/ 3> and the M — N number of line-solitons for 
{/ <SC 0. In particular, we give detailed description of the soliton solutions associated with Gr(2,4), 
which play a fundamental role in the study of multi-soliton solutions. We then consider a physical 
application of some of those solutions related to the Mach reflection discussed by J. Miles in 1977. 
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1. Introduction 

The purpose of this paper is to give a survey of our recent works on a classification theory of 
sohton solutions of the Kadomtsev-Petviashvili (KP) equation [H [T71 [31 El [7] . We also explain the 
detailed structure of some of the soliton solutions obtained in this classification theory, and discuss 
a physical application of a subclass of those solutions, which is related to the resonant interactions 
of the solitary waves in shallow water. Most of the materials in this paper are based on a series 
of lectures given by one of the authors (YK) at Chinese Academy of Science in Beijing after the 
international conference "Nonlinear Waves - Theory and Applications", on June 9-12, 2008. It is a 
pleasure to thank Professors Qing-Ping Liu, Xing-Biao Hu and Ke Wu for the arrangement of the 
lectures, and Mr. Kai Tian for making lecture notes. 

In 1970, Kadomtsev and Petviashvili [TS] proposed a two-dimensional dispersive wave equation 
to study the stability of one soliton solution of the KdV equation under the influence of weak 
transverse perturbations. This equation is now referred to as the KP equation, and considered to 
be a prototype of the integrable nonlinear dispersive wave equations in two-dimension. The KP 
equation is a completely integrable system with remarkably rich mathematical structure which is 
well documented in several monographs (see for examples, [H [H [H [T31 [501 [Ml UHl [2S])- These 
mathematical structures include the existence of A^-soliton solutions, the Lax formulation of the 
inverse scattering transform, and the existence of an infinite dimensional symmetries. However the 
real breakthrough in the KP theory occurred in 1981 by Sato [SI] who realized that the solution of 
the KP equation is given by a GL(oo)-orbit on an infinite dimensional Grassmann manifold (Sato 
universal Grassmannian) , and the KP equation is just a Pliicker relation on the Grassmannian. The 
present paper deals with a finite dimensional version of the Sato theory, and in particular, we are 
interested in the solutions which are real and non-singular in the entire xy-plane where they are 
localized along certain rays. 

It is quite important to recognize that the resonant interaction plays a fundamental role in multi- 
dimensional waves. The original description of the soliton interaction of the KP equation was 
based on a two-soliton solution found in the Hirota form, which has the shape of "X" describing 
an intersection of two lines with some angle and a shift of lines at the intersection (the phase 
shift). In 1977, Miles [25] pointed out that this two-soliton solution (referred to as the "0-type" 
solution, where "O" stands for original) becomes singular if the angle of the intersection is smaller 
than certain critical value. Since the KP equation was proposed to describe quasi-two dimensional 
waves, it should admit solutions that provide good approximations for two soliton interaction with 
smaller angle. Thus, it seems strange that one does not have a reasonable solution in the parameter 
regimes where the KP equation is supposed to give a better approximation. Miles also found that at 
the critical angle the two line-solitons of the 0-type solution interact resonantly, and a third wave 
(soliton) is created to make a "Y-shape" solution. Indeed, it turns out that such Y-shape resonant 
solutions are exact solutions of the KP equation (see also [26 ). 

After the discovery of the resonant phenomena in the KP equation, several numerical and exper- 
imental studies were performed to investigate resonant interactions in other real two dimensional 
equations such as the ion-acoustic and shallow water wave equations under the Boussinesq approxi- 
mation (see for examples [TBI [19l [27l [lOl [23] ) . However, after these activities, no significant progress 
has been made in the study of the solution space or real applications of the KP equation. It would 
appear that the general perception was that there were not many new and significant results left to 
be uncovered in the soliton solutions of the KP theory. For the last 5 years, we have been working 
on the classification problem of the soliton solutions of the KP equation [H [TT] [31 [El [^ , and our 
studies have revealed a large variety of soliton solutions which were totally overlooked in the past. 
In this paper, we give a brief survey of our research and explain an application of our new soliton 
solutions to describe the resonant interaction of solitons in shallow water. 

The paper is organized as follows: 
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In Section |2] we describe an interesting connection between the KP equation and the Burgers 



equation (Theorem 2.11, and then construct some exact solutions based on this connection. We 
show that the resonant interaction in the KP equation is a natural consequence of the mathematical 
structure of the KP equation, that is, the Burgers equation gives a symmetry of the KP equation in 
the y-direction. The Wronskian determinant structure of the solutions is derived from an extension 



of the Burgers hierarchy (Theorem 2.5 see also [T5] ). 

In Section |3] we present a brief introduction of the Sato theory, and give a direct construction 
of the Wronskian structure of the r-function (see also |3Ulllj). The KP hierarchy is the set of the 
symmetries of the KP equation, and the geometric structure of the solutions of the KP hierarchy is 
given by the r-function. We explain the important role that the r-function plays in the integrability 
of the KP equation. 

In Section|4] we give an elementary introduction of the real Grassmann manifold Gr{N, M), the set 
of iV-dimensional subspaces of M.^ . This is necessary background information for the classification 
problem of the soliton solutions of the KP equation. In particular, we explain the basic structure of 
the r-function in terms of the Schubert decomposition of the Grassmannian Gt{N, M). 

In Section 5] we present the main theorem of the classification problem of the soliton solution 
(Theorem 5.4). This theorem states that the r-function identified as a point on Gr(A^, M) generates 
a solution that has asymptotically M — N line-solitons for i/ ^ and N line-solitons for y 3> 
0. Moreover, these soliton solutions can be parametrized by the chord diagrams associated with 
certain permutation called derangements of the symmetry group Sm- This type of solutions is called 
(Af — N, A^)-soliton solution. We give a new proof of the Theorem, some of which did not appear 
in the previous papers. Specifically in Corollary |5.1| we show how the conserved densities of the 
KP equation lead directly to establish the correspondence between the line-soliton solutions and the 
derangements of Sm- 

In Section|6] we study the detailed structure of (2, 2)-solitons defined on Gr(2, 4) as a fundamental 
class of soliton solutions of the KP equation. We present the detailed description of the solutions 
based on the 2x4 matrix which marks the point of Gr(2, 4). The solution is completely determined 
if one prescribes a specific data for this matrix. We then show that for some classes of (2, 2)-soliton 
solutions, the matrix data is equivalent to asymptotic data such as the soliton locations and phase 
shifts, while for some other solution classes, one also need to specify internal data given by a local 
structure. 

In Section [7] we discuss an application of our new solutions to the Mach reflection in shallow water 
with a rigid wall, based on the results of Section [6] The resonant interaction among the incident 
wave, the reflection wave and the Mach stem is well described by our soliton solution, if we ignore 
the effect of the boundary layer at the wall. We also present a direct numerical simulation of the KP 
equation with a V-shape initial wave, which represents the reflection problem of the incident wave 
on an inclined wall (see also [12l [29j [33] ) . 

2. The KP equation 
The KP equation is the following partial differential equation in 2 + 1 dimensions, 

(2.1) {-Aut + 6uUx + Ua:xx)x + 3(3Uyy = 0. 

When (3=1, this equation is referred to as the KPII equation, while the equation with /3 = —1 is 
called the KP I equation. Since we consider only the KPII equation in this paper, we simply refer 
to the KPII as the KP equation throughout the text. 

Note when u is independent of y, the well-known Kortweg-de Vries (KdV) equation is obtained 



from (2.11 



-Aut + 6uUx + Uxxx = , 

which is why the KP equation is regarded as a (2 -I- l)-dimensional extension of the KdV equation. 
Physically, the KP equation was introduced to study the stability of the KdV soliton under the 
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influence of weak transverse perturbations So the KP equation should give a better approxi- 
mation of the original physical system if the system depends weakly on the transverse direction, i.e. 
the waves are almost parallel to the y-axis. 
Let w — w{x, y, t) be a function defined by 

u = 2wx , 



then the KP equation (2.1 ) takes the form 

{-AWxt + l2WxWxx + Wxxxx)x + SWxyy = . 

Integrating the above equation once with respect to x and setting the integration constant to zero 
(i.e. w is assumed to be constant as ~> oo), give the potential form of the KP equation, namely, 

(2.2) {—iWt + + Wxxx)x + SWyy ^ 0. 



We will henceforth refer (2.2 1 as the pKP equation for the function w which plays an important role 
in our description of the soliton solutions. 

We next review an interesting connection between the KP equation and the Burgers hierarchy, first 
established systematically in |T3]. This connection is the basis of our description of the fundamental 
structure of the soliton solutions discussed in this paper. 

2.1. The Burgers hierarchy and the KP equation. Let us first briefly review the Cole-Hopf 
transformation linking the heat (linear) and Burgers hierarchies: 

Proposition 2.1. /// is a solution of the set of linear equations, i.e. the heat hierarchy, 

(2.3) dtj^d'^f for n = l,2,..., 

then w :— dxQn. f) satisfies the following set of nonlinear equations, i.e. the Burgers hierarchy, 

(2.4) dt^w = dxidx + w^-^w for n = l,2, .... 
Proof. First, by induction, we prove that 

d:f=[idx + wr-'w]f:^Pnf. 

When n = 1, the above formula is just the definition of the function w. Now, suppose d^^^f — 
[{dx + wr-^w]f := F„_i/. Then, 

= dx{Pn-lf) = dx{Pn-l)f + Pn-iwf - [{dx + w)Pn-l]f = Pnf ■ 

Therefore, dt,w = a^J^.ln/) = dx{dtjf-^) = dx{d^^ff-') - 9.(9. + w^-^w. □ 

Note here that the Burgers hierarchy is well-defined in the sense that all the members commute, 
because the equations for / in (|2.3| are commute trivially. The Burgers equation corresponds to 



n = 2 in (2.4 1, i.e. the first nontrivial member of the hierarchy (2.4 1, and the members for n — 2,3 
are given by 

C dt^w = dx{wx + w"^), 



(2.5) 



dt^W ^ dx{Wxx+iwWx+W^) 



Note that the even members of the hierarchy are dissipative, while the odd ones are dispersive. 
This feature is a key for the resonance phenomena in the soliton solutions of the KP equation (see 
Example [O I . 

Now we show the connection between the KP equation and the Burgers hierarchy (2.4): 



Proposition 2.2. Suppose w(x,t2,t3) is a common solution of {2.5) with t2 = y and 3 — t, then 
w = w{x, y, t) solves the pKP equation ^2.^. 
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Proof. Expressing Wyy, Wxt in terms of only the x-partial derivatives of w, we have 

Wyy = d^{Wa;xx + '^wl + Aww^a: + 4w^W2:), 

Wxt = dxiwxxx + Swl + SwWxx + Sw'^Wx)- 

Then ehniinating the common terms including (wWxx + w'^Wx), i.e. calculate 3wyy — Awxt, we obtain 
the pKP equation □ 

Thus the KP equation contains the Burgers hierarchy, and the next theorem combines the results 
of Propositions |2. l] and pT2l which gives a linearization of the KP equation: 



Theorem 2.1. // /(x, t2, is, . . .) satisfies the set of linear equations, dt^f — f , n = 1, 2, . . . with 
t2 '■= y,t3 := t, then w — dx(l'a-f) satisfies the pKP equation (2.2). 

This implies that a solution of the KP equation obtained from the Burgers hierarchy (equivalently 
the set of linear equations) has a dissipative behavior in the y-direction. In particular, one should 
note that a confluence of shocks |Mj in the Burgers equation corresponds to a fusion of solitons as 
shown by Example |2.4| of the next section. 

Remark 2.2. One should also note that the infinitely many commuting symmetries of the Burgers 
hierarchy (2.4 1 induces higher flows of the pKP hierarchy. In other words, proceeding in the same 
way as in Proposition 2.2 one could also obtain an equation for w{x,t2^ ■ . ■) in the first n variables 
x,t2,t3, . . . ,tn, which form the nth member of the pKP hierarchy. The substitution u = 2wx then 
gives the corresponding member of the KP hierarchy. Thus, the KP equation (2.1) corresponds to 
the n = 3 fiow of the KP hierarchy, and the higher fiows generate the infinitely many symmetries 
of the KP equations. The KP hierarchy plays an important role when we discuss the multi-soliton 
solutions as we show in this paper. 

2.2. Some exact solutions. We now consider some solutions of the KP equation obtained from 
the linear system (2.3 1. The general solution for this linear system with t2 ~ y and = t admits 
the integral representation (Ehrenpreis principle) 

f{x,y,t)^Je'^+'''y+''''p{k)dk, 
c 

with an appropriate measure p{k) dk and a proper contour C in the complex plane. A particularly 
simple finite dimensional solution is recovered by choosing 

M 

p{k) dk = a.j6{k — kj) dk 
with arbitrary real constants kj and aj for j = 1, 2, . . . , M, and the contour C = M. Then 

M 

(2.6) f{x,y,t)=Y.ajE,{x,y,t) E,ix,y,t) := e''^-'y'*\ 

with 0j := kjX + k'jy + k^t. Here we assume aj to be positive, so that / is positive definite, and also 
assume the ordering in the fc-parameters, 

(2.7) ki <k2<--- < kM. 
Now let us give some explicit examples: 

Example 2.3. Consider the case with Af = 2, i.e. / = aiEi + a2i?2. Since w = dx{ln f), we can 
choose ai = 1 without altering the solution. So we take, with oi = a, 

f = Ei + aE2 = 2Vae5(«i+»2) cosh ^(6*1 - 6*2 + 612), 
where 612 = — In a. Then we have 

w = 292 In/ = i(fci - k2)hech^{9i -62 + ^12) , 
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Figure 2.1. A One-soliton solution. The [1,2] is the label of this line-soliton, and 
each E(j) for j — 1,2 indicates the dominant exponential term Ej in this region. 
The amplitude A[i,2] find the angle ^'[1^2] s-rc given by ^[1,2] — ^(fe ~ ^1)^ ^'^d 
tan^'[i 2] = fci + fc2. K[i 2] = {K^i 2]'-^[i 2]) shows the wave-vector, and the slope 
is given by tan^[i_2] — Kj^ 2]/ ^[i 2]- Throughout the paper, the graphs of soliton 
solutions illustrate contour lines of the function u{x,y,t). 



which is called the one-soliton solution. The solution u(x, y, t) is localized in the a;y-plane along the 
line 61 — 62 + 9i2 = 0, whose location is determined by the constant a. In particular, at t = the 
peak of the soliton determines the line, 

x + {ki + k2)y = - — ^—r- In o, 
fc2 - fci 

so that choosing a = 1 the line of the peak crosses the origin. Because of this, we often refer the 
one-soliton as a line-soliton (in a local sense in general) . 

We also note that in terms of the function w{x, y, t), we have the following asymptotic values for 
x — > ±00 with the ordering fci < ^2, 

fx kiEi + ^2-^2 f fc2 as X ^ +00, 

w — — = ^ < , 

/ Ei + E2 { ki a.s X —00. 

and the solution u = 2wj. 0, exponentially as |a;| — > 00. The change in the asymptotic value of the 
potential function w for the one-soliton solution is due to fact that the exponential Ei dominates 
over E2 SiS X —^ —00, whereas the exponential E2 dominates over Ei as x 00. Consequently, the 
solution u is localized along the phase transition line: 9i — 62 ^constant, where both exponential 
terms are in balance. This motivates labeling the one-soliton solution as the [1, 2]-soliton which also 
represents a permutation ti" = (2 1) of the index set {1,2}, i.e. exchanging the values ki and k2 by 
crossing the line-soliton. 

It is also convenient to introduce the following notations to describe a Zme-soliton physically as 
a traveling wave: We say that a line-soliton is of [i, jj-type (or simply [i, j]-soliton), if the soliton 
solution u has the form (locally), 



(2.8) 



u = sech^ (K[,j] • x + U^.^^^t + ^l^^^ 
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The amphtude ^[ij], the wave-vector K[jj] and the frequency ^[ij] are defined by 

1 

1 



1 



1 



2"(% — ki) {l,ki + kj) , 



The direction of the wave- vector K[i j] = (iiTj^ , if ^ ) is measured in the counterclockwise from the 
X-axis (see Figure 2.1 ), and it is given by 



tan* 



Note here that the ^'[ij] also gives the angle of the line from the y-axis in the counterclockwise, 
and — I < ^'[i.j] < f . For each soliton solution of (2.81, the wave vector K[j jj and the frequency 
f2[ij] satisfy the soliton-dispersion relation, 

(2-9) 4J7[,,,-jif[X,1-4(Xp,,])' + 3(if[l,])'. 

The soliton velocity V[i jj defined by ■ = given by 



n 



K, 



K, 



■1,3] I 



ko k'i "4" k- 

' ^(_i, _(fc, + fc^.))• 



l + (fci + kj)^ 



Note in particular that {kf + kikj + kj) > 0, and this implies that the x-component of the velocity is 
always negative. That is, a line-soliton propagates in the negative x-direction. On the other hand, 
any small perturbation propagates in the positive x-direction. This can be seen from the dispersion 
relation of the KP equation for a linear wave (/> = exp(ik • x — icut) with wave-vector k = (fc^,, ky) 
and frequency uj, 

1 3 S ky 

4 4: k^ 

from which the group velocity of the wave is given by 
V = = 



du) duj 
Qkf dku 



3 



Thus, it is reasonable to expect from a first-order perturbation theory that the soliton separates from 
small radiations asymptotically, similar to the case of the KdV equation (see for example [T| [25]). 

We also remark that the formula (2.8 1 for one-soliton solution of the KP equation can be extended 
to the one-soliton solution of the KP hierarchy by including the higher times (see Remark 2.2), 

(2.10) u{h,t2,t3,...)=Asech'' JC.ntn + Q° 

with some constant 0". Here the amplitude and the infinite dimensional wave- vector are given by 

A = 2^^-' ~ 

/C„ = ^(fc^-fef), for n=l,2,.... 
This can be easily seen from the structure of /-function, i.e. 

f^E, + aEj, with E, = exp ( £ k^t^ 

\n=l 
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Although the higher times t„ for n > 3 do not have a direct physical meaning, those parameters are 
related to the existence of the multi-soliton solutions through the symmetries of the KP equation as 
mentioned in Remark 2.2. 

Example 2.4. Now we consider the case with M = 3. We again take ai — 1, so that we have 

with some positive constants a and b. As in the previous example, it is also possible here to 
determine the dominant exponentials and analyze the structure of the solution in the xy-plane. Let 
us consider the function / along the line x = —cy with c = tan^* where 5* is the angle measured 



counterclockwise from the y-axis (see Figure 2.11. Then along x = —cy, we have the exponential 
function Ej ~ exp[rij{c)y + kjt] with 

(2.11) T]j{c)^kj{kj-c). 



It is then seen from Figure pT2| that for y ^ and a fixed t, the exponential term Ei dominates when 
c is large positive (^' « ^, i.e. x — oo). Decreasing the value of c (rotating the line clockwise), 
the dominant term changes to E^. Thus we have 



w = dx In f 



ki as X — oo, 
ks as a; — > oo. 



The transition of the dominant exponentials Ei — > E^ is characterized by the condition 771 = ?73, 
which corresponds the direction parameter value c = tan^'ji 3] = ki + k^. In the neighborhood of 
this line, the function / can be approximated as 

f^Ei + bEs, 

which implies that there exists a [1, 3]-soliton for y ^ 0. The constant b can be used to choose a 
specific location of this soliton. 

Next consider the case of y <C 0. The dominant exponential corresponds to the least value of 
rjj for any given value of c. For large positive c (^' ~ ^, i.e. x 00), E^ is the dominant term. 
Decreasing the value of c (rotating the line x = —cy clockwise), the dominant term changes to E2 
when k2 + ks > c > ki + k2, and Ei becomes dominant for c < ki + k2- Hence, we have for y <C 

{ki as X ^ —00, 
k2 for -{ki + k2)y < X < -{k2 + k3)y, 
fca as X ^ 00. 

In the neighborhood of the line x + (fci + k2)y =constant, 

f^Ei+aE2, 
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Figure 2.3. A (2, l)-soliton solution. Each E{j) with j — 1,2 or 3 indicates the 
dominant exponential term Ej in that region. The boundaries of any two adjacent 
regions give the line-solitons indicating the transition of the dominant terms Ej. 
The fc-parameters are the same as those in Figure |2.2| and the line-solitons are 



determined from the intersection points of the ?7j(c)'s in Figure 2.2 Here a = b 



(i.e. T = Ei + E2 + E3) so that the three soUtons meet at the origin a,t t — 0. 



which corresponds to a [1, 2]-soliton and its location is fixed by the constant a. The solutions also 
consists of a [2, 3]-soliton in the neighborhood of the line x+{k2 + k3)y =constant, and whose location 
is determined by the locations of other line-solitons. Therefore, we need only two parameters a, b 
(besides the fc-parameters) to specify the solution uniquely. These free parameters may be considered 
as the asymptotic data for the solution u. The shape of solution generated by / = i?i -I- aE2 + bE^ 
with a = b = \ (i.e. at i = three line-solitons meet at the origin) is illustrated via the contour 
plot in Figure p3| In this Figure, one can see that the line-soliton in ?/ ^ labeled by [1,3], is 
localized along the phase transition line 9i — 9^ (equivalently, 771 = 773) with direction parameter 
c = ki + ks; two other line-solitons in y <C labeled by [1,2] and [2,3] are localized respectively, 
along the phase transition lines with c ~ fci + fc2 and c = ^2 + ^3. This solution represents a resonant 
solution of three line-solitons. In terms of the function w, which is a solution of the Burgers equation 
in the y-direction, this corresponds to a confluence of two shocks (see p. 110 in |31])- The resonant 
condition among those three line-solitons is given by 

K[i,3] = K[i^2] + K[2,3], ^[1,3] = ^[1,2] + ^[2,3] I 

which are identically satisfied with H-nj] = -^{kj — ki, — kf) and ^[i.j] = ^{kj — kf). The resonant 
condition may be symbolically written as 

[1,3] = [1,2] + [2,3]. 

One can also represent this line-soliton solution by a permutation of three indices: {1,2,3} which 
is illustrated by a (linear) chord diagram shown below. Here, the upper chord represents the [1,3]- 




soliton in y ^ and the lower two chords represent [1,2] and [2, 3]-solitons in y ^ 0. Following the 
arrows in the chord diagram, one recovers the permutation. 



1 2 3 
3 1 2 



or simply tt — (312). 
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Figure 2.4. The time evolution of a (3, l)-soliton solution and the corresponding 
chord diagram. The upper chord represents the [1, 4]-soliton, and the lower ones 
represent [1,2]-, [2,3]- and [3, 4]-solitons. The chord diagram shows tt — (4123). 



In general, each line-soliton solution of the KP equation can be parametrized by a unique permutation 
corresponding to a chord diagram. We will discuss this issue in Section 5. 

The results described in the previous examples can be easily extended to the general case where 
/ has arbitrary number of exponential terms (see also \21\ |4]). 

Proposition 2.3. If f = aiEi + 02-^2 + • • • + a^EM with aj > for j = 1,2, .. . ,M, then the 
solution u consists 0/ M — 1 line-solitons for ?/ <C and one line-soliton for y 3> 0. 

Such solutions are referred to as the (M — 1, l)-soliton solutions. Note that the line-soliton for 
J/ ^ is labeled by [1,M], whereas the other line-solitons in y <C are labeled by [fc, fc -I- 1] for 
k — 1, 2, . . . , Af — 1, counterclockwise from the negative to the positive x-axis, i.e. increasing \E' from 
— "I to |. As in the previous examples one can set ai = 1 without any loss of generality, then the 
remaining M — 1 parameters a2, . . . ,aM determine the locations of the M line-solitons. Also note 
that the a:?/-plane is divided into M sectors for the asymptotic region with ~\- 3> 0, and the 
boundaries of those sectors are given by the asymptotic line-solitons. This feature is common even for 



the general case. Figure 2.4 illustrates the case for a (3, l)-soliton solution with / = Ei+E2+E^+E4^. 
The chord diagram for this solution represents the permutation vr = (4123) of the set {1,2,3,4}. 



2.3. Extension of the Burgers hierarchy. Theorem 2.1 shows the relation between the KP 
equation and the Burgers hierarchy for the function w defined in terms of the solution / of the linear 



system (2.3 1 via the Cole-Hopf transformation, 
(2.12) 



w 



dx In /, 

where /^^-^ — dxf ■ Furthermore, Proposition 



/ 



(1) 



wf, 



2.3 



illustrates the {M — 1, l)-soliton solutions obtained 
from the special choices for the functions / and w. In order to construct more general type of soliton 



solutions of the KP equation, we consider an extension of the Cole-Hopf transformation (2.121 to 
the case where the function / satisfies an iV-th order linear equation, 

(2.13) /W ^ + • • • + WN^if^'^ + w^f, 

where f^^^ := d^f. The coefficient functions {wi,W2, ■ ■ ■ ,wn) can be constructed from N linearly 



independent solutions of (2.131, and their time evolutions with respect to n > 1 are determined 
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by the compatibility conditions of ( 2.13[ ) with the hnear system ( |2.3[ ). For example, when N = 2, the 
compatibility condition dy{fxx) = with y = t2 gives dy{wif^^^ + W2f) = (9^(wi/(^' + ^2/)- 

Then by equating the coefficients of /(^' and / in the last expression, leads to the following equations 
for {wi,W2), 



dyWi 

dyW2 



2W2Wl.a; 



Wl.xx 
W2.XX 



■ 2W2,: 



Notice that if W2 = 0, this system is reduced to the Burgers equation. In general, the variables 
(wi, W2, . . . ,wn) satisfy an N -component Burgers equation in the y-variable while the higher flows 
with respect to the "times" t„ for n = 3,4, .. . are the symmetries of this coupled equation, thus 
forming an iV-component Burgers hierarchy (this will be reformulated in terms of the Sato theory 
in the next section, also see jl 3j ) . 

We regard (2.131 as an A'^-th order ordinary differential equation, and consider a fundamental set 
of N solutions denoted by {/i, /2, . . . , /at}, which satisfy 

+ + , = 1,2,. ..,7V. 

The above form a linear algebraic system which can be solved for the coefficients {wi, W2, ■ ■ ■ , wpf} 
by applying Cramer's formula. In particular we have 



1 

Wr(/i,...,/^) 



= |:lnWr(A, 



fi 

/2 

In 
Jn) 



(1) 



(1) 



JN 



f{N-2) 

AN-2) 
J2 



AN-2) 
JN 



(N) 



(N) 



AN) 
JN 



where Wr(/i, . . . , /jv) is the Wronskian determinant of {/i, . . . , /at}. For = 1, this is just the 
Cole-Hopf transformation (2.121 where wi = w solves the pKP equation. In what follows, we show 
that this is also true for the general case. That is, the function wi — dx lnWr(/i, . . . , /jv) solves the 
pKP equation, and u := 2wi^x = 29^ In Wr(/i, . . . , /jy) solves the KP equation. The Wr(/i, . . . , f^) 
is an example of the r-function for the KP equation, which plays a very important role of the 



KP theory (see Section 3.3). We now state the well-known theorem that leads to the Wronskian 



formulation of the multi-soliton solution for the KP equation. 

Theorem 2.5. Let the r-function be given by the Wronskian determinant, 



(2.14) 



fi 

/2 



/i 



(1) 



(1) 



/i 
f: 



{N-l) 



{N-1) 



N 



JN 



AN-1) 
JN 



where {/i, /2, . . . , /n} o,i~e linearly independent functions satisfying dyfj = /j^'' and dtfj — fj^ with 
/■"^ :— d!^fj for j — I, . . . , N . Then u{x, y, t) ~ 28^ lnT(x, y, t) satisfies the KP equation {2.1). 



We prove the theorem by employing an identity called the Pliicker relation (see for example |14j 
satisfied by the maximal minors of an N x M matrix 



^ := 



/2 /: 



(1) 



(Af-l) 



N 



/•(I) 

Jn 



AM-l) 

Jn 
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with N < M . A maximal minor ip{li, I2, ■ ■ ■ , In) is the determinant of the submatrix formed by N 

columns = {ft'\ft'\ ' ' ' : Jn^'Y of ^ where ie{h,l2,-- - , In}, i.e. 

ifHuh, . . . , ^at) = det [^i, ,<Pi„..., . 

Lemma 2.1. T/ie sei 0/ maximal minors (p[i, j] :— (p{l, 2, ■ ■ ■ , N — 2, N — 2 + i, N — 2 + j) of the 
matrix <P defined above, with 1 < i < J < 4, satisfy the identity, 



(2.15) 



(^[1, 2] ^[3, 4] - (^[1, 3] ^[2, 4] + (^[1, 4] ^[2, 3] = 0. 



Proof. Consider the following 2N x 2N determinant expressed in terms of the columns <l>i of the 
matrix (p, 

'<Pl ... <Pn-1 ^1 ••• 'Z'JV+I <Z'Ar+2" 

... ^1 ... f?Ar_2 ^AT <?Ar+i ^N+2 

This determinant is obviously zero. Then, Laplace expansion of the first determinant in terms of 
N X N minors yields the desired identity. □ 



det 



The three-term identity (2.151 is one of the Pliicker relations which play an important role in 



revealing the geometric structure underlying the KP equation. This will be discussed in Section |4] 
where we will introduce the general form of the Pliicker relations. We now prove Theorem |2.5| (the 
following proof is the same as that in [Mj): 

Proof. Substituting the relation u = 29^(lnr) into the KP equation (2.1 1, integrating twice with 



respect to x and setting the integrations constants to zero, yields the following equation in r. 



(2.16) 



4(rrj.t - T^Tt) - 3t|^ - rr^ 



3Ty - iTTyy = 0. 



The r-function is given by the maximal minor, r = ip{l, 2, . . . , N) =: ip[l, 2], using the notation of 



Lemma 2.1 Suitable combinations of the derivatives of r also correspond to other maximal minors 
of the matrix (p. In particular, we have 

■^[1,3] =T^, ip[l,i] = l{T^a: + Ty), (p[2,3] = ^{t^.^ - Ty) , 



Tt), ip[3, 4] = ]^(Tj,2::rj; + 3Tyy ~ 4r^t) 



Substituting the above expressions into the Pliicker relation (2.151 in Lemma 2.1 we have 
=^[1,2] ip[3, 4] - ^[1,3] ip[2, 4] + ^[1,4] (p[2, 3] 



= -1^(^(4^:. 



.)) 



which is equivalent to ( |2.16[ ). □ 

Theorem |2 . 5| provides a direct linearization scheme for the KP equation in the sense that a large 



class of solutions whose r-function is given by the Wronskian determinant (2.14) can be constructed 
from the solutions of the linear equation (2.13) together with the evolution equations fy = fxx and 
ft — fxxx- Note that since f(x,y,t) admits a rather general integral representation (see Section 
2.2), this linearization scheme gives rise to different classes of solution besides the soliton solutions, 
by proper choice of measure and contour of integration in the fc-plane. It is a remarkable fact that 
all solutions constructed in this way identically satisfy the KP equation by virtue of the Pliicker 



relation (2.15) which is an algebraic identity. However, it should be emphasized that the significance 



of Theorem 2.5 is much deeper than that is presented here. In the original scope of Sato theory, the 



differential operator in (2.13) is in fact a reduction (A^-truncation) of a pseudo-differential operator 
of infinite order (see Section 3) all of whose coefficients are determined via a single holomorphic 
function namely, the Sato r-function. The Wronskian determinant (2.14) is a representation of the 



Sato r-function in the A^-truncated case. It is in this (infinite-dimensional) setting of the full Sato 



theory that the KP equation expressed in terms of the Sato r-function as in (2.16) can itself be 



interpreted as the Pliicker relation (2.15) in a suitable sense (see also Remark 4.4) 



Remark 2.6. The equation (2.16) is often rewritten in the Hirota bi-linear form 
P{Dx,Dy,Dt)T ■ T := {ADxDt - - 3Dl)r . r = 0, 
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where with z = x, y or t is the usual Hirota derivative defined by 



9e" 



[f{z + e)g{z-e)] 



'=0 



Note here that the soliton dispersion relation ( |2.9[ ) can be obtained from the Hirota bi-linear form, 
i.e. P{2K^^jj,2K^^ , ^^liij] ) = 0. Thus the dispersion relation has a direct connection to the Pliicker 
relation (2.15). 

In his famous book [14 , Hirota expresses the minors j] with the elegant diagrams called the 
Maya diagrams introduced by Sato (which is also equivalent to the Young diagrams). Explicit proofs 
of the solution formulas for other integrable systems written in the Hirota bi-linear forms are also 
given in terms of the Maya diagrams. 

Example 2.7. Consider the case with N = 2 and Af = 3: The linearly independent functions /i 
and /2 are expressed as, 

3 

fi= J2 O-ij^j 1 * = 1,2, 



with Ej{x, y, t) = exp(kjX + k'jy + fc|i) for j 
explicitly given by 



1,2,3. In this case the r-function in (2.141 can be 



A n 

f2 f, 



(1) 



(1) 



Oil ai2 

021 ^22 ^23 





kiEi\ 


)p 


k2E2 




k^E^j 



ail ai2 

021 122 



i?(l,2) 



On 

021 



ai3 
a23 



E{1,3) 



0.12 ai3 

122 0-23 



where we have used the Binet-Cauchy theorem and E{i, j) ~ Wi{Ei, Ej) 
investigate a concrete situation where {ki,k2,k3) ~ (— |,— i, |), and 



E{2,3), 

{kj — ki)EiEj. Let us 



On 012 

021 ^22 



ai3 

023 



1 
1 



with positive constants o and b. The r-function is then given by 

E{1,2)+ aE{l,S) + bE{2, 3). 

In order to carry out the asymptotic analysis in this case one needs to consider the sum of two T]j{c), 
i.e. r]i j — rji + rjj for I < i < j < 3- This can still be done using Figure 2.2 but a more effective 
way is described in Section 5 (see the graph of ry(fc, c) in Figure 5.1 and equations (5.4|, (5.5l). 
For y 3> 0, the transitions of the dominant exponentials are given by following scheme: 

E{1,2) ^ Eil/S) E{2,3), 

as c varies from large positive (i.e. x — > —oo) to large negative values (i.e. x — > cxi). The boundary 
between the regions with the dominant exponentials E{1,2) and i?(l,3) defines the [2, 3]-soliton 
solution since here the r-function can be approximated as 



E{1,2) + aE{l,3) = (fc2 - ki)Ei { E2 

1 



ki 



2+03 



'^^^ ) cosh i( 6*2 - 



so that we have 



where ^23 is given by 



2(fc2-fci)£;ie2 
u = 29| Inr w i(fc2 - A;3)2sech 
ki 



723 



In 



ki 



In a 



I.e. 



k2 — ki 
+ O23), 

; + ^23) , 
-fci 



E. 
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Figure 2.5. A (1, 2)-soliton solution. The E{i,j) on each region indicates the 
dominant exponential term. The /c-parameters are the same as those in Figure |2.2[ 
The parameters in the v4-matrix are chosen as a — ^ and & = 1, so that three 
line-solitons meet at the origin at t ~ 0. 



Note here that since the logarithm of the coefficient of cosh- function in r above is linear in x, it 
does not contribute to the solution u. The constant a can be used to determine the location of this 
soliton. A similar computation as above near the transition boundary of the dominant exponentials 
E{1, 3) and E{2, 3) yields 

r w 2{h - fci)a£;3e5('i+'=-'^=) cosh -02 + ^12), 
2 1 , 



^^{ki - k2fscch^{9l - 02 + 012) 



where 0i2 is given by 



. ks-ki b . k2~ki _ 
In ; in — I.e. = — e 

fca - fc2 a k^ - fc2 



This gives [1, 2]-soliton and its location is determined by the constant b. 
For J/ <C 0, there is only one transition, namely 

i?(2,3)^i?(l,2), 

as c varies from large positive value (i.e. x — > 00) to large negative value (i.e. x — *■ —00). In this 
case, a [1, 3]-soliton is formed for j/ <C at the boundary of the dominant exponentials E{2, 3) and 



i?(l, 2). The contour plot of the line-soliton solution is shown in Figure 2.5 Notice that this figure 



can be obtained from Figure 2.3 by changing {x,y) {—x, —y). This solution can be represented 



by the chord diagram corresponding to the permutation tt — (231) shown below. Note that this 




diagram is the n-rotation of the chord diagram in Example 2.4 whose permutation tt = (312) is the 
inverse of tt = (231). 

3. Sato theory of the KP hierarchy 

In this section we review briefly some basic facts from the Sato theory of the KP hierarchy (see 
also [HI [25 )• fi'^st review the Lax formalism for an infinite set of evolution equations whose 
compatibility conditions give rise to the fiows corresponding to the KP hierarchy. The KP equation 
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is the first nontrivial member of this hierarchy associated with the time variables ti := x, t2 '■= y 
and := t. The main purpose is to highlight the basic framework of integrability underlying the 
KP theory. Then we show that the results of the previous section are direct consequences of the 
Sato theory. In particular, we emphasize the importance of the r-function and explain the central 
role of the r-function in the KP hierarchy. 

3.1. Lax formulation of the KP equation. Let L be a first order pseudo-differential symbol, 

L=^d + ui + U2d^^ + usd^'^ H 

where the coefficients Ui = Ui(t) depend on infinitely many variables t = (ii, t2, t^, ■ ■ ■). The symbol 
d := dx is a differential operator whereas is a formal integration, satisfying dd^^ = d~^d = 1. 
The operation of 9" with G Z is given by the generalized Leibnitz rule 

]>0 V J / 

For example, we have Of — fx + fd and d^^f = fd^^ — fxd^^ + fxxd^^ • (the latter expression 

follows from the usual integration by parts formula). Here we define the weights for the functions 
Ui and as 

wt(Mi) — i, wt(d'^) — v. 
So the L has a homogeneous weight one. We also remark that the term ui in L can be eliminated 
by a gauge transformation g such that ui = —g^^dx{g), i.e. 

L — > g-^Lg = d + ua^"^ + u^d-^ + ■■■ . 
(This can be extended so that we can eliminate all Uj with an appropriate gauge operator g = W, 



which will be discussed in Section 3.2 ) We will henceforth consider L without the Ui term, 

(3.1) L^d+U2d-' + U3d-^ + --- . 

The Lax form of the KP hierarchy is defined by the infinite set of nonlinear equations 

(3.2) at„L=[S„,L] with S„-(L")>o n=l,2,..., 

where (L")>o represents the polynomial part of L" in d, i.e. Bn is a differential operator of order n, 
and [Bn, L] := BnL — LBn is the usual commutator of operators. Since [S„, L] = [L" — (L")<o, L] = 
[L, (L"')<o], each side of the equation (3.2 1 is a pseudo-differential operators of order < —1. Here 
(L")<o is the negative part of in d, and note [d, (i")<o] has no polynomial part. Thus for n > 1, 
each equation in (3.2 1 is a consistent but infinite system of coupled (1 + l)-evolution equations in tn 
and X, for the variables {ui : i = 2,3, ■ ■ ■}. The case n = 1 yields the equations dt^^Ui — dxUi, so we 
identify <i with x. The infinite system (|3.2| is compatible, as prescribed by the following theorem: 



Theorem 3.1. The differential operators i?„ = (L")>o satisfy 
(3.3) dt^Bn - dt^Brn + [B„, S„,] = 0. 

Consequently, the flows defined by (3.^ commute i.e., for any n,m > 1, 



Proof. It follows from (3.2) that dt^{L") = [Bjn,L"]. Then using the decomposition L" 
Bn + (L")<o, we have 

9t„(i") - dt^Ln = - [Bn,L"'] 

= [Bm,Bn] ~ [(L"')<0, (i")<o]. 



The differential part (> 0) of the above equation gives (3.3 1. 

To prove that the flows commute we compute using (3.2) once again, and obtain 

dt„Sdt„L) - dtSdt,M = [9t„S„,L] + [Bn,dt„L] - [9t„B^,L] - [B,n,dtM 

= [(^t„,Bn — dt,^B,n, L] + [Bn, [Bni,L]] — [i?,„ , [i?„, i]] 



16 



SARBARISH CHAKRAVARTYt AND YUJI KODAMA* 



Applying the Jacobi identity for commutators, the right hand side of the above equations becomes 
[dt„^Bn — dt^B„i + [Bn, B„i],L], which vanishes due to (3.3 1. □ 

Equations (3.3) are called the Zakharov-Shabat equations for the KP hierarchy. Note that for a 
given pair (n,77i) with n > m, (3.3) gives a closed system of n— 1 equations for U2,M3, • • • in the 
variables tm,tn and x. For example, consider the case with n = 3 and m = 2, i.e. B2 = { L^)> q 



+ 2u2 and -B3 = (L'^)>o = 9'^ + 3u2 + 3(u2,a; + W3), then the Zakharov-Shabat equation (3.3) for 
B2 and B^ gives the system 

f dt^U2 = U2.XX + 2U3^X 

I 29(3 U2 = 3{u2,x + U3)t2 - {U2,xx - 3u3,a; + 3^2)2; 

After setting <2 = y,t3 — t, and eliminating U3 from the system, u — 2u2 satisfies the KP equation 
(pll). 



We also remark that the KP hierarchy (3.2) is given by the compatibility of the linear system 
(3.4) 



L(l> 



= Br, 6 



1,2, 



where fc € C, the spectral parameter, and (f>{t;k) is the eigenfunction of the KP hierarchy. The 
compatibility among the second set of equations gives the Zakharov-Shabat equations. 



3.2. The dressing transformation. As we mentioned in the previous section, the Lax operator 
L can be gauge-transformed into the trivial operator d, i.e. 

(3.5) L — >d^ W-^LW, 
where the operator of the gauge transformation is defined by 

(3.6) W = I - wid-^ - W2d-'^ ~ w^d^^ + ■ ■ ■ . 



The coefficient functions Wi are related to the coefficients Uj of L via the relation LW 



Wd, and 



which yields 

U2 



Wj^x 



U3 = W2,x + WiWi^x , "4 = W3^x + {wiW2)x 
+ Fj + i{wi,W2,...,Wj^l), ••• , 



wj ^ + wlwi^x ■ 



Thus the Wj can be 



where i^j+i are differential polynomials of weight j -I- 1 (note wt{wj) — j) 
considered as the primary variables whose x-derivatives determine the KP variables as illustrated by 
the KP and pKP equations mentioned in Section 2.1 The evolutions of the Wj with respect to the 
time variables t„ can be prescribed in a consistent fashion by requiring that the dressing operator 
W satisfy the following system of equations. 



(3.7) 



dt^W = BnW-Wd'' 



for n — 1,2, 



where i?„ is now given by Bn = {Wd^W~^)>Q. Notice that this expression for Bn as a differential 
operator is a consequence of the equations, [(9t^ VF)M^^^]>o = 0. Equation (3.7) is sometimes 
referred to as the Sato equation. 

The following Theorem asserts that the KP hierarchy is generated by the dressing (or conjugation) 
of the trivial commutation relation [dt^ — 9", 9] = by the operator W. Because of this result, the 
(inverse) gauge transformation, d ^ L, is called the dressing transformation for the KP hierarchy. 



Theorem 3.2. If the oper ator W satisfies the Sato equation (3.7), then the operator L = WdW~^ 
satisfies the Lax equation (3.2) for the KP hierarchy, and the operators Bn — {Wd^W~^)>Q = 
(L")>o satisfies the Zakharov-Shabat equations (3.3). 
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Proof. First, a direct calculation using L — WdW ^ and i" — Wd^^W ^ = B„ 
that the operator 



(i")<o shows 



where the last equality is due to (3.7). Then equations ( |3.2[ ) and (3.3 1 follow from the commutator 
relations 



= W[dt„ - 5", d]W-^ - dt„L - [Bn, L] , 

= W[dt^ - 5", dt^ - d^W-^ = - B„, dt, 



B„ 



□ 



It follows from Theorem 3.2 that the flows defined by the Sato equation (3.7) are commutative 



i.e. they satisfy the compatibility condition dt,^{dt^W) — dt„{dt„^W). Indeed, the compatibility 
condition for ([^ is equivalent to [dt,^ + (L")<o"9tJ+ (L'")<o] = 0. Using L" = B„ + (L")<o, the 
commutator term on the left hand side can be decomposed as [dt^ — Bn, 9t„ — Bm] + [dt^ — Bn, i'"] — 
[^tm ~ ^m, i"] , which vanish due to Theorem (3.2 1. Finally we note that the KP linear system (3.4) 
is obtained by the dressing action: (j) — 
linear system 

(3.8) ^ ^ ''^'^ 



W(j)Q where the vacuum eigenfunction (j)Q satisfies the hare 



1,2, 



This equation together with the Sato equation (3.7) form the basic ingredients of the dressing 



procedure. We will use the bare eigenfunction in the normalized form. 



(3.9) 



with 



(t;fc) = ^ k'Hr, 



3.3. The T-function. The r-function (2.14) introduced in Section 2.3 plays a fundamental role in 



Sato theory for the KP hierarchy. In this subsection we demonstrate explicitly how the r-function is 



related to the dressing operator W which satisfies the Sato equation (3.7). We restrict our discussion 



in this paper to a finite truncation of the infinite order pseudo-differential operator W for simplicity 
only, while capturing the flavor of the general theory. Let us then consider the dressing operator 

W^l- wi(9^i - W2d-^ wnQ-^ , 

and deflne the differential operator, 

Wn Wd"" = 9~ - wid""-^ - W2d''-^ WN . 

Since W satisfles the Sato equation, it is immediate that Wn also satisfles 



dt„ 



= BnWn - WnO"^ Bn = {Wd''W-') 



>o 



The following Proposition establishes the compatibility conditions leading to the extended Burgers 
hierarchy introduced in Section |2] 

Proposition 3.1. The N-th order differential equation Wjsif = is invariant under any flow of the 



linear heat hierarchy (2.3), {dt^f = d"f : n = 1,2, . . .}. 



Proof. It suflices to show that dt^ (WOv/) = 0. Then the desired result follows from the uniqueness 
of the differential equation: 

dt^WNf) = dtAWN)f+WNdtJ 

= {BnWN~WNd")f + WNdtJ 
- Bn{WNf) + WN{dtJ^d':f)^0. 
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Proposition |3.1| provides the compatible system considered in Section |2.3 
(3.10) 



r WNf = 0, 

\dtj = d^f, 



71= 1,2, . 



Furthermore, a set {fj : j = 1, 2, . . . , N} of hnearly independent solutions of W^f = 0, i.e. 

can be employed to explicitly construct the coefficient functions Wi of the dressing operator W as 
follows 



(3.11) 



Wi 



(-1) 



i+l 



fl 
/2 



/i 



(AT-i-l) .(AT-i + l) 



N 



AN- 



AN- 
JN 



fl 
J2 



-1) AN-i+1) 
J N 



/i 



(AT) 



AN) 
J N 



where r Wr(/i, . . . , /at). Note that since the i„-dependence of the Wi is given via the evolution 
equations dt„fj — d^fj ior j — 1,2, . . . , N, one also has an explicit solution of the Sato equation 

The formula (3.11 1 for the coefficients Wi can be exploited to obtain an elegant expression for the 
KP eigenfunction <f) via the r-function (see for example [51 l24jV 

Lemma 3.1. The eigenfunction cf) of the linear system (g.^P can be expressed as 



(t;fe)=lil_g:^^o(t;A:) 



r(t) 



where 4>o = e^**''^) with Oit; /c) = E k'X,, 



(t-[fc-i]) -.^{t,--, t, 



2k'- 



3fc3' 



We provide a simple proof of this well-known formula to preserve the self-contained style of the 
present article (see also [5]). 

Proof. First note that using (3.111, the eigenfunction can be written as 



= W^0o = 1 - 



fc2 



fl 



(1) 



fl 



JN 



(N) 



fN 

k-N j^-N+1 ... 



/i 



AN) 
JN 



Using the elementary column operations, the determinant in the numerator of the above expression 
can be re-written as 

(-1)^ 



(^"-'^"■') ...... 



From the integral representation of the functions fi. 



/,(t) = /e'^(*■^V^(A)dA for z = 1, 2, . . . , iV , 
c 
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each matrix element in this determinant is given by 



c 

where we have used ln(l — f) = — S^""- This completes the proof. □ 

n=l " 

One should note that the expression for (j) in Lemma 3.1 also holds in the general case for the 
full untruncated version of the operator W. Expanding this formula with respect to k, we have an 
explicit formula for Wi's in terms of the r-function, i.e. 

Wi = Si{-d)T, 

T 

where d :— {dtn ^^tj, \dt.^, . . .) and S'„(z)'s are the complete homogeneous symmetric functions (or 
sometimes referred to as the elementary Schur polynomials) defined by 



(3.12) exp^ fc".2„ = ^ ^„(z)fc", z = (zi,Z2,...)- 

\n=l / n=0 

The explicit expressions of these polynomials are given by 

S'o(z) = 1, S'i(z) = Zi, 5'2(z) = Z2+y, ^3 (z) = Z3 + Z1Z2 + y , ... , 

I1+2I2A hn/„=n 'l''2! ' hi- 

Note that for an A^-truncated operator W (i.e. w„ = if n > N) the r-function satisfies the 
constraints 

S'„(-a)T = for n>N. 

Lemma |3.1| will be now used to derive a set of first integrals of the KP equation that will prove 
to be useful in our classification of the line-soliton solutions in Section 5. The integrability of the 
KP equation may be demonstrated by the existence of an infinite number of conservation laws in 
the form, 

'^tnhj = dxgj.n , 

for some conserved densities hj and the corresponding conserved fluxes These functions are 

differential polynomials of Wj's in the Lax operator L, and they can be found as follows: Differenti- 
ating the quantity (f)~^dx<t' with respect to t„ and using the evolution equation dt„(t> = Bn(j>, we first 
derive the conservation law, 

(3.13) 9t„ {cl)-'dxcf>) = dx [r^Br^cf) . 



Next we invert (3.1 1 using the generalized Leibnitz rule to express the differential operator d in dx4) 
in terms of L as 

d = L- V2L~^ - v^L^'^ , 

where vj's are the differential polynomials of Mi's. Then the conserved density 4'~^9x(f> can be written 
as an infinite series after using L^cf) — k"(j), n G Z, 

a,(/.= fc---p , 
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Each function Vj is a conserved density of the KP hierarchy; the first few are given by 



Note that these are nonlocal in the KP solution, u = 2u2 (except ^2). However, they can be expressed 
simply as certain differential polynomials of the variable wi = ^^^(Inr). Namely, we have 



^dx(j) = d^ncl}^ [6'(t;fc) + lnr(t - [fc-^]) -lnT(t 

1 



exp 



n=l 



1 



1 



InT(t) = A;+ T^d^Sni-d) Inr, 

n=l 



Vn+1 = -dxSn{,-d) Inr = -S'„(-9) wi 



which leads to 
(3.14) 

For example, we have 

V2 dl In T, i;3 = \dx {dt^ - dl) In r, Vi = \dx {dt^ - ^d^dt^ 

Vn+i = (<9t„ + (h.o.d.)) Inr, • • • . 

where (h.o.d) indicates the terms including higher powers of the derivatives. Moreover, if the solu- 
tions Ui's of the KP hierarchy decrease rapidly to zero as |a;| — > 00, one can define the integrals C„ 

by 



i5^)ln. 



(3.15) 



Cn ■■= J Vn+i{x, ...)dx 



for 



1,2, 



In particular, if the r-function gives one line-soliton of [i,j]-type, i.e. T(t) — Ei(t) + aEj(t) with 
Ei{t) — exp(6'(t; ki)), then from (3.141 the integrals take 

1 



(3.16) 



-Sni-d) Inr 



Notice here that if ki < kj, t k Ei ^ (,e(t;ki) for a; < and r « aEj = ae^(*'''j) for x » 0. In 



general, if there are N line-solitons of [i;, j;]-type for I — 1, . . . , iV, then we have 

k?) , 



N I 

1=1 n 



3l 



(see Corollary 5.1 below for the proof). This expression is similar to the A''-soliton solutions of the 
KdV equation (see for example [33]), but for the KP equation, the C„'s are also independent of y. 
This fact will be used in Proposition |5.1| below. 



Remark 3.3. An alternative set of the conserved densities can be found by observing the following 
tautological equations (see p. 99 in [S]), 

dtrr,{dxdt„ Inr) = dx{dt^dt„ Inr). 

That is, the conserved densities obtained from those equations are given by 

Vn+i ■■= -dxdn Inr. 
n 

^ CO 00 

However one can show that the integrals C„ = / Vn+i dx = J Vn+i = C„ for all n. 



In Section |5] we will discuss the classification problem of the solutions obtained from the r- 
function with finite dimensional solutions of the /-equation, and will show that the classification 
is completely characterized by the asymptotic behavior of the r-function. Before getting into the 
classification problem, in the next Section, we give a brief review on the real Grassmann manifold 
for the background information necessary for the problem. 
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4. Introduction to the Grassmannian Gr{N, M) 

The main purpose of this Section is to explain a mathematical background of soliton solutions of 
the KP equation which is provided by the real Grassmannian Gr(Af, M), the set of iV-dimensional 
subspaces of M*^ This is a finite dimensional version of the universal Grassmannian introduced by 
Sato in 1981 [31] • 

4.1. Grassmannian Gy{N^ M). Let {/i, /2, • • ■ , /at} be a set of N linearly independent vectors in 
M^^. Then the set forms a basis for V , 

^ = SpanR{/i,/2,...,/jv}, and Gr(iV, M) = : F C K^} . 
Let us begin with some simple examples: 

Example 4.1. The Grassmannian Gr(l,Af): Suppose M — 2, then Gr(l,2) is the set of all lines 
through the origin in the plane. Thus, any two non-zero vectors in are on the same line if they 
are proportional, and they determine the same point of Gr(l,2). Since every point (a;i,X2) ^ (0,0) 
of lies on some line I through the origin, we can represent I, which is also a point on Gr(l,2), as 

Gr(l, 2) 9 / = (xi : X2) {c{xi,X2) | c ^ 0, {x^x^) € \ {(0, 0)}} . 

The line {xi : X2) is called the homogeneous coordinate for Gr(l,2) which has then the structure of 
a one-dimensional manifold. This manifold admits the decomposition, 

Gr(l,2) = {(1 : a) : a e M} U {(0 : 1)} = MU {00}, 

which is the projective line MP^, and it can be identified with the unit circle (i.e. a e M in the 
decomposition parametrizes the slope of each line, and the y-axis is identified as (0 : 1)). 

The set of all lines passing through the origin of is represented by the Grassmannian Gr(l, M) 
which can be identified with the M — 1 dimensional projective space MP*^^^. Homogeneous coor- 
dinates and the decomposition of Gr(l,Af) can be described in the same way as in the previous 
example. In this case, one has the well-known Schubert decomposition given by 

M 

Gr(l, M)=\J F.J with = {(0 : . . . : : 1 : ai : . . . : au-j) ■ Ofc € K} , 

where each cell Fj has the structure of a (M — j)-dimensional vector space; the last one being 
Fm = (0 : . . . : : 1) which represents a point (the 0-dimensional cell) in Gr(l, M). 

We now consider a representation for the general case of Gr(A^, M) which plays a key role in the 
classification of the line-soliton solutions of the KP equation as we show in the next section. With 
respect to a basis {Ei,E2, ■ ■ ■ ,Em} of M*^, the basis vectors fi of the iV-dimensional subspace V 
can be expressed as 

M 

(4.1) f,= Y.a^jEj, for i = l,2,...,N, 

i=i 

where A :— {uij) is a, N x M matrix. Note that the rank of A is iV since the {/i, /2, ■ • ■ j /at} are 
linearly independent. A different choice of basis for V will induce row operations on the original 
matrix: A BA with B E GL{N). Consequently, any matrix A can be canonically chosen in the 
reduced row echelon form (RREF) by taking appropriate B G GL(iV). This canonical form of A has 
a distinguished set of pivot columns labeled by X = {zi, 12, ... , in} with 1 < ii < i2 < ■ ■ ■ < in l£ M 
such that the submatrix Aj formed by the column set 2 is the identity matrix Ipf. For any given 
A, its unique RREF provides a coordinate for a point of Gr(A^, M). Moreover, the set of all such A 
matrices in RREF, and with a fixed pivot set X forms an affine coordinate chart or cell Wx which 
(similar to Gr(l, M)) gives the Schubert decomposition of the Grassmannian 

(4.2) GiiN,M)= y Wj. 

IC{1,...,A/} 
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Each cell is called a Schubert cell. For example, ifX — {1,2,..., N}, then the Schubert cell Wx 
contains all A matrices whose RREF is given by 



/ 1 
1 



V 



* 
* 



where the N{M — N) entries of the right-hand block are arbitrary. This particular Scubert cell Wx 
is often refered to as the top cell of the Grassmannian Gr(A^, M), and has the maximum number of 
free parameters marked by *. It follows from this that the dimension of Gr(A^, M) is N{M — N). 
A straightforward counting of the maximum number of free parameters for an A in RREF with a 
given pivot set X = {ii, . . . , shows that the dimension of the cell Wx is given by 



dim Wx = N{M - N) + ^N{N + 1) - (n + + 



Example 4.2. For Gr(2,4), the Schubert decomposition (4.2 1 is given by 



Gr(2,4) 



u 



W, 



There are six cells W^ijy with dim W{i.j} 
(a) W{i,2} = 

(c) W^{1,4} = 
(e)H^(2,4>^{(g ; 



l<i<j<4 

7 — (i + j), which are given by 

1 i 
( 



(b) W^{i,3} = 

(d)H^(2.>^{(;; I 

(f) W^{3,4} = 





1 



1 
1 





4.2. The Pliicker embedding and the r-function. We now introduce a coordinate system called 
the Pliicker coordinates for the Grassmannian Gr(iV, AI) which is related directly to the Wronskian 
form of the r-function (2.141 in Theorem |2.5| Given any basis {/i, . . . ,/Ar} of an iV-dimensional 



subspace V C M*^, the exterior product space /\^ is a one-dimensional subspace /\^ V C /\^ 

/i A /2 A • • • A /n- If {51, ■ . ■ , gw} is another basis for V with 

(A, ^ - - — - 



spanned by the wedge product lyj 
B e GL(iV) so that (51, . . . , .gAr) 
an element of P{/\^ M.^'^), the projectivization of /\ 



fN)B, then Vg = det{B)h'f. Consequently, /\ is 
The Pliicker embedding is a map from 
M.^ . In terms of the A-matrix in 
(4.1), the explicit form of this map is given by 



Gr(7V,M) to P(A^M^), that sends V C R^' to /\^ V C /\ 



(4.3) 



/i A /2 A • • • A /, 



N 



l<ii<---<iN<M 



where ^(ii,i2, ■ ■ • ,^Af) are the maximal minors of A formed by its columns Ai, i e {ii,i2, • ■ ■ , iw}- 



Recall that the maximal minors were introduced in Section 1.3 (following Theorem 2.5). The minors 
in the set {^(ii, • • • , lAf) : 1 < ii < • • • < «Ar < M} are defined up to a scale factor, and are called 
the Pliicker coordinates for the Grassmannian Gr(A^, M). These coordinates are not independent, 
but satisfy the identities 



(4.4) 



N+l 

L (-l)'~^C("i7 • • • ,aAf-i, A)C(/3i, • • • , • • • ,Pn+i) = 0, 



called the Pliicker relations which hold because the exterior form i^f — fi A ■ ■ ■ A /n is totally 
decomposable, i.e., v Avf = 0, y v €z V . The Pliicker relations can be derived using elementary linear 
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algebra from the Laplace expansion of the following 2N x 2N determinant formed by the columns 
Aj of the matrix A, similar to Lemma 12. 11 



Aai 





Af-l 





Ab. 



0. 



The Pliicker coordinates, modulo the Pliicker relations, give the correct dimension of Gv{N^ M) 



which is typically less than the dimension of P{[\ 



Example 4.3. For Gr(2,4), the Pliicker coordinates are given by the maximal minors, 

e(l,2), e(l,3), e(l,4), C(2,3), e(2,4), ^(3,4). 
Taking ai — 1, {Pi^ Pi) — (2,3,4) in ( |4.4| gives the only Pliicker relation in this case, 

2)^(3, 4) - ^(1, 3)^(2, 4) + 4)^(2, 3) - 0, 



which is the same as (2.151. Since dim(/\ M*) = 6, and the projectivization gives dim(P(A M*)) = 
6 — 1 = 5. Then with one Pliicker relation, the dimension of Gr(2,4) turns out be 4, which is 
consistent with the dimension of the top cell W{i^2} as shown in Example 3.2 (case (a)). 



We next show that the r-function defined by the Wronskian determinant in ( 2.14 ) can be identified 
as a point on the Grassmannian Gr(Af, M). First note that the set of functions. 



exp kpr 



J 



1,2, 



with distinct real parameters kj, gives a linearly independent set because Wr{Ei, E2, ■ ■ ■ ,Em) 7^ 
0. Then an iV-dimensional subspace V C M.^ is given by Spanjj{/i, /2, . . . , /n} with N linearly 
independent solutions {/i, /2, . . . , /at} of fy — f^x and ft = fxxx as in Theorem 2.5 The coefficients 



of Ej in (4.1 1 form a, N x M matrix A of rank N, whose row-space is isomorphic to V so that A 
represents a point of Gr{N, M). One can then see that the r-fu nction r = Wr(/i, /2, . . . , /at) gives 
the Pliicker embedding Gr(iV, M) — > P{/\^ M.^^) defined in (4.3 1: Indeed expanding the Wronskian 
determinant (2.14) by Binet-Cauchy formula, we have 

(4.5) 



Wr(/i,/2,...,/jv) 



J2 C(«l,»2, ■ • ■ ,«jv)£^(il,»2, ■ 
l<il<---<ilv<A/ 



■ ,«Af) 



where S.{ii,i2, ■ ■ ■ ,iN) are the Pliicker coordinates given by the maximal minors of the A-matrix, 
and E{ii,i2, ■ . . ,«Ar) = Wr{Ei-^, Ei^, . . . ,Eij^). Note here that E{ii,i2, ■ ■ ■ ,ijv) can be identified as 



Ei, A E,, A 



■ Eij^ , and forms a basis for /\ 



NmM 



N 



Spanjj {E{ii,i2, ■ ■ ■ Jn) ■ ^ < ii < i2 < ■ ■ ■ < iN < M} = f\] 



Note here that the sum fc^^ 



kij^ should be distinct for distinct sets {11,12, . . . , Jat} in 



order for the {E{ii,i2, . . . , iat)} to be independent. 

Remark 4.4. We just mention here that an infinite dimensional Grassmann manifold appears when 
we consider the case where the functions {/i, . . . ,/Ar} are analytic. In this case, one can expand 
fi{t) in terms of the elementary Schur polynomials defined in (3.121, 



Pi{k)dk^ ai,„S'„_i(t), 



c 



where 9(t; k) = Y ^"tn, and ai.„ are given by the moment integrals. 



ai,n ^ J k"- pi{k) dk 
c 



for 



l<i<N, 
Kn. 
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Thus the A-matrix has the size N x oo, and each A-matrix represents a point on an infinite dimen- 
sional Grassmann manifold Gr(A^, oo). The r- function then has the well-known expansion form with 
the Schur polynomials, 

(4-6) T(t) = J2 C(ji,i2, • ■ • ,iAf)5'{ji,,2,...,,„}(t), 

l<ii<i2<-<ijv 

where S[i-^^,,,^i^j(t) is defined by 

5'{ii^j2:---,ijv}(t) = Wr(S'ji_i,S'i2_i, . . . ,S'i„„i)(t). 
The index set {11,12, ■ ■ ■ ,*Ar} is commonly expressed by a Young diagram (see for example [24'). 



Then the main result of the Sato theory is that the function defined by (4.6 ) is the r- function, which 
represents a point on the Grassmann manifold |31j . 

In the following section, we will show that the line-soliton solutions arise from the r-function given 



by (4.5). The Grassmannian formulation of this r-function provides a geometric characterization of 
the solutions by identifying the different solution classes with distinct Grassmann cells. Of particular 
physical interests are those solutions which are regular in x, y and t. Such solutions are realized if 
the associated r-function has non-negative Pliicker coordinates, that is, all ^{ii,i2, ■ ■ ■ ,iN) > in 



(4.5). The matrices A having this property are called totally non-negative (TNN) matrices; they 



parametrize the totally non-negative Grassmannian Gr+(iV, Af) C Gr(iV, M) (see also [5U]1. 

5. Classification of soliton solutions 
In this section, we present a classification scheme of the line-soliton solutions based on the asymp- 



totic behavior of the r-function (4.5). It turns out that the classification problem is closely related 
to the characterization of the totally nonnegative Grassmannian cells of Gr(iV, Af ) recently studied 
by Postnikov et al (see for example [30| [35]). 

5.1. Asymptotic line-solitons. Let us recall that the r-function is given by the sum of exponential 
terms, 

(5.1) T{t) ^ J2 i{mi,. . . ,mN)E{mi,. . . ,mN){t), 

l<mi<---<mjv<Af 

where ^(mi, . . . , mjv) = Wr(£;,„,, ... ,£;„„) with E„,{t) = e«(*;'='") = exp( g fc;;,t„). We then 
impose the following conditions: 

(i) The /c-parameters are ordered as ki < k2 < ■ ■ ■ < kM , and the sums ki + kj are all distinct 
(this implies that all [i, j]-solitons have different slopes). 

(ii) The ^-matrix is TNN, that is, all its maximal minors ^("^1; "^2i ■ ■ ■ are non-negative 
(i.e. the r-function is positive definite for all x,y,t, hence the solution u — 2c)^(lnr) is 
non-singular) . 

The asymptotic spatial structure of the solution u(x,y,t) is determined from the consideration 
of which exponential term Eimi, . . . ,712^) in the r-function dominates in different regions of the 
xy-plane for large \y\. For example, if only one exponential term E(mi, . . . ,misi) in the r-function 
is dominant in a certain region, then the solution u = 29^(lnr) remains exponentially small at all 
points in the interior of any given dominant region, but is localized at the boundaries of two distinct 



regions where a balance exists between two dominant exponentials in the r-function (5.1). 

Furthermore, it is also possible to deduce the following result from the asymptotic analysis of the 
r-function: 

Proposition 5.1. The dominant exponentials of the r-function in adjacent regions of the xy-plane 
are of the form E{i,m2, ■ ■ ■ ,rnpf) and E(j,m2, . . . ,miq). That is, the exponentials contain N — 
1 common phases and differ by only one phase. 

We prove the Proposition by applying the following Lemma. 
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Lemma 5.1. Let {s„(xi, X2, ■ ■ ■ , xj^) : n — 1, . . . , N} be the set of N power sums in N real variables 
defined by 

s„ := x'/ + x^' H hx^, n = 1,2, . . . ,7V. 

Suppose the set of N equations 

Sn{xi,X2,- . ■ ,Xn) ^ Cn, n^l,...,N, 

admit solutions (xi, . . . , xjv) for appropriate values of the constants c„, then the solution set is unique 
up to the permutation of N elements {xi, . . . , xn}- 

Proof. The proof follows from the Newton's identities relating the power sums {s„} and the 
elementary symmetric polynomials {cr„} which are defined as 

CT„ = J2 Xi^Xi^---Xi^ n=l,...,N . 

l<ii<---<in<N 

For n = 1,2, ... ,N, the Newton's identities are recursively given by 

Sn + aiSn-1 + . • • + a„-isi + nan = , aj, = (-l)''crfc . 

Then the solutions of s„ = c„, n — 1, . . . , N, which exist by assumption, correspond the set of N 
roots xi, . . . ,xn of the iV-th degree monic polynomial, 

x^ + aix'^-^ + a2X^-^ + • • • + aw = 0, 

where the coefficients {a„} are uniquely determined by {c„} using Newton's identity. Obviously, the 
set of roots {xi, . . . , x„) is unique up to permutations. □ 
Now we prove Proposition |5.1| 

Proof. Let E{ii, . . . , z^r) and E{ji, . . . ,jN) be the two dominant exponentials in adjacent regions, 
then near the boundary of the two regions, the r-function can be approximately given by 

T w £,iE{ii, . . . ,?Ar) + • ■ • Jn), 

where ^1 and ^2 £^re the corresponding minors of the A-matrix. The remaining terms in the r-function 
are exponentially small in comparison with those dominant terms. Then the solution u = 29^(lnT) 
locally has the one-soliton form as in (2.101, i.e. 

/ 00 

\n=l 

where A= ^iQ-P)'^ and JCn = ^(Q" - P"). Moreover, P and Q satisfy 

Q" _ P« ^ (fc« + . . . + fc^J - {kl + • • • + fcr„) n^l,2,.... 

Here it is assumed (without any loss of genericity) that kj-^ + • • • + kj^ > ki-^ + ■ ■ ■ + ki^ so that 
P < Q. Let us rewrite the first + 1 of the above equations as follows 

(5.2) P» + fc« + . . . + = Q« + fc^'; + . . . + fc^^ := c„ , 

where c„'s are values of each sum. Then {P, kj^, . . . , kj^ } and {Q, ki^, . . . , ki^ } are solutions to the 
system of equations x" + . . . x^j^^ — n — 1, . . . ,N + 1. Hence, from Lemma 



5.1 



these two sets 

must be the same up to permutation of their elements. Now recall that the sets {fci„} and {fcj„} 
are distinct and P ^ Q. Then the one-to-one correspondence between the sets {P,kj-^, . . . ,kj^} 
and {Q, ki-^, . . . , ki^^} implies that P must be one of {ki^^}, and Q must be one of {fcj„}. Then the 
remaining iV — 1 elements of the set {ki^} must be the same as the remaining N —1 elements of the 
set {fcj„}, proving the Proposition. □ 

As an immediate consequence of Proposition |5.1| the asymptotic behavior of the KP solution is 
given by 

(5.3) u{x,y,t) « i(fc, - h)hech^^{9j -0, + %) , 

in the neighborhood of the line x + (ki + kj)y— constant , which forms the boundary between the 
regions of dominant exponentials E(i, ni2, . . . , m^) and E{i, m2, . . • , rn^). Equation (5.3 1 defines 
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Figure 5.1. The left figure shows r]j{c) = kj{kj — c) for j = 1, . . . ,4. Each at 
the intersection point of 77^(0) — r]j(c) indicates the exchange i <~+ j. We assume that 
there is at most one intersection point for each c (genericity condition for kj 's) . The 
right one illustrates the 77 as the function of k with c — fc2 + ^3, that is, the 77 along 
the dotted line in the left figure passing through the intersection point 772(c) = 774(c). 
This figure shows the order 773 < 772 = 774 < 7^1. 



an asymptotic [7, j]-soliton as a result of those two dominant exponentials. Then the condition (i) 
below (5.1 ) implies that those asymptotic solitons have the different slopes, so that they all separate 
asymptotically for \y\ ^ 0. 

Next we discuss how to determine which exponential term in a given r-function is actually dom- 
inant in the direction of particular line x = ~cy for y — > ±00. The basic idea is the same as in 
Example 2.4 Along the line x — —cy, each exponential term E(mi, 777,2, ■ • ■ , wat) has the form, 

/ N 

Cexpij: 77(fc™„,c)y + C„W 



where C and O^^it) are constants for fixed t, and ri{k,c) is defined by (cf. (2.11 1), 
(5.4) 



77(A;, c) — k{k — c) 



and 



77j(c) = kj{kj -c). 



, c) for each c. First note 



Then for y (or ^ 0), we look for the dominant (or least) sum of 77(fc, 
that we have 

r]i{c) = i]j{c) if c = ki + kj, 

that is, in the direction of [i, j]-soliton, those terms are in balance. Since the fc-parameters are 
ordered as /ci < fc2 < • • • < k^J, we have the following dominance relation among the other 77^(0) 's 
along c — ki + kj, 



(5.5) 



Vj < Vm if 1 < 77T, < 7 or J < 771 < M , 

Vj > Vm if i < m < j. 



In order to find the dominant sum, the graph of ri{k,c) is particularly useful. Figure 5.1 illustrates 
the case with M = 4. We demonstrate how to identify the asymptotic line-solitons from a given 
T- function using Lemma 5.1 and the relations (5.5). 



Example 5.1. We consider the A-matrix given by 

A 



110 
11 



Then from (5.1 1 the r-function has the form, 

T^E{1,3)+ E{l,i) + E{2, 3) -I- E{2, 4). 



SOLITON SOLUTIONS OF THE KP EQUATION 



27 



Since Proposition |5.1| implies that the hne-sohtons are locahzed along the phase transition lines 
X + cy =constant with c = ki + kj , we look for dominant exponential terms in the t- function along 
those directions. 

For y 3> 0, the c-values along a line x = — decrease as we sweep clockwise from negative to 
positive a;-axis. So, we start with the largest value c = k^ + k^, and try to determine if the line 



[3,4] actually corresponds to a line-sohton. The relation (5.5 1 gives ?7i,7y2 > % = Vii along the line 
[3,4]. This implies that r/i + 772 is the most dominant combinations along this line. However, since 
the corresponding exponential term £'(1,2) is absent in the r-function, it turns out that i?(l,3) 
and £(1,4) are the two dominant exponentials in the r-function along the line [3,4] satisfying the 



condition of Lemma 5.1 Therefore, there is a [3, 4]-soliton. An almost identical argument implies 
that [1,2] is another line-soliton. But for example, along the line [1,3], the following relations hold: 
> Vi — Vst V2 < Vi = V3- This implies that r]i + r]4 and 7?3-|-t?4 are the two dominant combinations. 
But £(1,4) is the only dominant exponential term along the line [1,3], since £(3,4) is absent from 
the r-function. Hence the solution u « 0. In this fashion by checking along the lines for all 
possible (j,j) pairs, we conclude that [1,2] and [3,4] are the only two asymptotic line-solitons as 
2/>0. 

For y <C 0, we now look for the least sum rji + rjj which implies the dominant exponential 
E(i,j). We continue to sweep clockwise from positive towards negative x-axis with the c-values still 
decreasing. Again, we start with c— fcs -I- fc4, then we have the relation ?/i, 772 > ~ Vi- -But since 
£(3,4) is absent from the r-function, £(2,3) and £(2,4) provide the dominant balance leading to 
the line-soliton [3,4]. Considering c = fc2 -I- k^, we have rj^ < r)2 = rji, 774, which implies that £(2, 3) 
is the only dominant exponential in the r-function, hence m « along the line [2,4]. Proceeding 
similarly as in the case for y ^ leads to the conclusion that [1,2] and [3,4] are the only two 
asymptotic line-solitons for y <C 0. 

From the above example, it is important to notice the fact that whether or not the dominant 



(or least) combinations inferred from (5.5 1 are actually present in the given r-function ultimately 
decides if a line corresponds to a line-soliton. In turn, this depends on the coefficients of the 
exponential terms given by the maximal minors of the A-matrix. Therefore, it suffices to specify 
only the A-matrix rather than the whole r-function in order determine the asymptotic line-solitons. 
The next example illustrates this point. 

Example 5.2. Consider the 2x4 matrix, 

A - 
1 + + 

where some of the non-zero entries are indicated by their signs which ensure that all non-zero 
maximal minors of A are positive. In this case, there are five nonzero minors, 

e(l,2), e(l,3), e(l,4), C(2,4), e(3,4), 

and one missing ^(2, 3). 

Now look for the asymptotic solitons as y 0. First we can immediately see that [3, 4]-soliton 



is impossible since the dominant exponents along the line [3,4] from (5.5) is given by rji + 7/2 



Since C(l;2) 7^ 0, T{x,y,t) ~ £(1,2) implying that u ~ along [3,4]. For the same reason. 



[1,4] and [1, 2] -solitons are also impossible. Let us then check the [2, 4]-soliton. From (5.5l, 773 < 
??2 = ?74 < Vij and since £,{^,2) ^ 0, ^(1,4) 7^ 0, the r-function has the dominant phase balance 
r w £(1, 2) H- ^(1, 4)£(1, 4) along [2,4]. Therefore [2,4] corresponds to an asymptotic line-soliton as 
y 00. Moreover, the [1, 3]-soliton also exists for similar reasons. Thus, we have two asymptotic 
line-solitons for y ^ 0. 

We next look for the asymptotic solitons for y <C 0. It is easy to see that [1, 2] and [3, 4]-solitons are 



impossible. Then consider the [1, 3]-soliton. In this case, (5.5 1 implies that rj2 < rji = rj^ < 774. Thus 
the two exponents £(1,2) and £(2,3) are dominant for y <C 0. However, ^(2,3) = implies that 
this is impossible. So [1, 3]-soliton does not exist as y ^ 0. For similar reasons, [2, 4]-soliton is also 
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Figure 5.2. A (2, 2)-soliton solution. The left left figure is at i = -16 and the 
right one at t — 16. Each E{i,j) for 1 < i < j < 4 in the left figure shows the 
dominant exponential in the region, and the boundaries of those regions give the 
line-solitons. The parameters (fci, . . . , ^4) are chosen as (—1, — ^, 2). 



impossible. Now check [1 , 4]-soliton. In this case, we have 772 , V3 < Vi = Vi from (5.51 (see also Figure 



5.1 1. Since ^(2,3) ~ 0, the dominant exponent E{2,3) is not present in the r-function, but there 
does exist a dominant balance with i?(l,3) and -B(3,4). Thus there is a dominant phase transition 
along the line [1,4]. A similar argument applies for the transition line [2,3] which corresponds the 
other line-soliton as ?/ ^ 0. In Figure |5T2] we show the soliton solution given by this matrix. 



5.2. Characterization of the line-solitons. We saw in the Examples |5.1| and |5.2| that the A- 
matrix plays a key role in the identification of the asymptotic line-solitons since they determine if 
particular dominant phase combinations are in fact present in the given r-function. Therefore, in 
order to obtain a complete characterization of the asymptotic line-solitons it is necessary to consider 
the structure of the N x M coefficient matrix A. For the remainder of this article, we will consider 
the matrix A to be in RREF, and we will also assume that A is irreducible as defined below: 

Definition 5.3. An N x M matrix A is irreducible if each column of A contains at least one nonzero 
element, or each row contains at least one nonzero element other than the pivot once A is in RREF. 

The reason for this assumption is the following: if an iV x M matrix A is not irreducible, then the 
solution u = 29^(lnr) can be obtained from a r-function with a matrix A of smaller size than the 
original one (that is, the size of A is reducible). If a column of A is identically zero, then it is clear 
that we can re-express the functions ji in terms of a x (Af — 1) coefficient matrix A obtained from 
A by deleting its zero column. Or, suppose that an A'^ x Af matrix A in RREF has a row whose 



elements are all zero except for the pivot, then it can be deduced from (4.5) that the corresponding 
r-function gives the same KP solution u which can be obtained from a another r-function associated 
with a (A^ — 1) X (M — 1) matrix A. For example, consider the matrix. 



A = 



1 ai2 ai3 
1 



which yields /i = £'1 aY2E2 + a^zE^ and /a £'4. Then r = Wr(/i,/2) = E^(k^h - h.x)- 
Factoring out the exponential £4, we find that both r and f = — fi^x give the same KP solution 
u. Furthermore, note that r = / = (i?i, £2, Ez)A^ , where A is the 1x3 matrix, 

A = (an, ai2, 013) = {k^ - fci, k^^- /c2ai2, ^4 - ^3013) ■ 

We now present a classification scheme of the line-soliton solutions by identifying the asymptotic 
line-solitons as y ^ ±00. We denote a line-soliton solution by (A^_, A^+)-soliton whose asymptotic 
form consists of A_ line-solitons as y — > — cx) and A^+ line-solitons for y — s- 00 in the xy-plane 



as shown in Figure 5.3 The following Proposition provides a general result characterizing the 
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Figure 5.3. (Ar_, iV+)-soliton solution. The asymptotic linc-solitons are de- 
noted by their index pairs [e„,j„] and [imj5m]- The sets {ei, 62, . . . , e^v} and 
{ffij 525 • ■ • 5 5A/-Af} indicate pivot and non-pivot indices, respectively. Here N_ = 
M — N and Af+ = N for the r-function on Gr{N, AI), and E{-, ...,•) represents the 
dominant exponential in that region. 



asymptotic line-solitons of the (Af_, Af_|_)-soliton solutions. 

Proposition 5.2. Let {ei, 62, . . . , cat} be the pivot indices, and let {51, 52, • ■ ■ , ffM-w} be the non- 
pivot indices of an irreducible and totally non-negative N x M matrix A. Then the soliton solution 



generated by the r-function in (5.1) with the matrix A has 

(a) N asymptotic line-solitons as y 00, each defined uniquely by the line [e„, j„] for some jn, 
and 

(b) M — N asymptotic line-solitons as y —00, each defined uniquely by the line [im,gm] for 
some im- 

Proof. Let [i,j] denote an asymptotic line-soliton for y 3> 0, and let it correspond to the dominant 



exponentials E(i,m2, . . . jTOjv) and E(j,m2, . . . ,m]^) in (5.11. Hence in particular, the coefficient 



minor ^(i,rn2, ■ ■ ■ tTUn) 7^ 0. Assume that Ai is not a pivot column. Then Ai must be spanned 

n 

by the pivot columns of A to its left since A is in RREF. That is, Ai — J2 CrA^^, where e„ < i. 

r=l 

Therefore, 

n 

^(i, m2, . . . , rriN) = J2 Cr£.{er, "^2, ■ • ■ , "^^v) 7^ . 

r=l 

Suppose Cs is the coefficient of the first nonzero term in the above sum, then the minor ^(ej, m2, . . . , mj^) 
0, implying that the corresponding exponential E{es, m2, . . . , tojv) is present in the r-function. Fur- 



thermore, from (5.5 1 rje^ > rji along so that E{es,m2, • . . , rnj^) > E{i, m2, . . . , tojv)- But this is 
impossible as E(i, TO2, . . . , wjv) is the dominant exponential in the r-function along [i,j]. Therefore, 
Ai must be a pivot column, proving part (a) of the Theorem. 

Now suppose [z, j] is an asymptotic line-soliton for y <C and that it corresponds to the dominant 
exponentials E{i, TO2, . . . jUIn) and E{j, 7712, . . . , tun) in ( |5.1[ ). Then the nonzero minor ^(i, 7712, . . . , ttiat) 
which implies that the column set B — {Ai,Am2, ■ ■ ■ ,A„ij^} forms a basis for M.^ . In particular, 

N 

B spans the column Aj, i.e., Aj — cAi -\- ^ Cr-Amr- Then, 

r=2 

^(i,m2, . . . ,771r_l,j, 771^+1, . . . ,7777v) = Cr^(7,77l2, . . . , 77ljv) , r = 2, . . . , iV . 
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If Aj is a pivot column, then there is at least one index s such that > j. Then it follows from (5.5 1 
that rjm^ > rjj along the line and for y <^ 0, we have E{i,m2, ■ ■ ■ ,ms-i, j,ms+i, . . . , tun) 3> 

E{i, 7712, . . . , rriiq). But this contradicts with E{i^ m2, . . . , tojv) being dominant. 

To prove uniqueness of each index pair, first consider ?/ > 0, and assume that there are two 
line-solitons [en,jn] and [e„, j,'J with a given pivot index e„, and where j„ > j'n- These line-solitons 
are localized along the lines x + {k^^^ + kj^)y = constant and x + (fcg^ + kji^)y = constant with 
^e„ + fcj„ > fce„ + kji^. Notice from Proposition 5.1 that the phase 9e„ is replaced by the phase 



in the exponential term which is dominant immediately to the right of the [e„, j„]-soliton. On 
the other hand, 0e„ must be present in the exponential term which is dominant immediately to 
the left of the [e„, j'^]-soliton. Consequently, there must be an intermediate line-soliton [em,e„] 
for some pivot index < e„, localized along the line x + (fce^ + ke„)y = constant, which must 
satisfy k^^ + fce„ > k^^ + ky . But the latter inequality is impossible due to ordering of the indices 
Sm < e„ < j'^ and the ordering kr < kg for r < s, giving the required contradiction. Hence for a 
given pivot index e„ there is a unique line-soliton [e„, as y — * oo. The uniqueness proof for the 
[im, 5m]-solitons for m — 1, . . . , M ~ N, and as y — > — oo, is similar. □ 



The unique index pairings in Proposition 5.2 have a combinatorial interpretation. Let [M] := 



{1, 2, . . . , M} be the integer set and recall that {ei, . . . , e^} U {gi, . . . , qm-n} is a disjoint partition 
of [M]. Define the pairing map tt : [M] [M] according to parts (a) and (b) of Proposition 5.2 as 
follows: 



(5.6) 




= j„, n = 1,2, . . . ,iV, 

= im, m= 1,2,..., M-iV, 

where e„ and gm are respectively, the pivot and non-pivot indices of the A-matrix. Then as a 
consequence of Proposition |5.2| we have the following. 

Corollary 5.1. The map tt : [M] [M] is a bijection. That is, n G Sm, where Sm is the group of 
permutations for the index set [M] . 

Proof. We prove this by finding an explicit relation among those parameters {fci„,%„} and 
{fce^jfcg^}: For this purpose, we consider the integrals C„ in (3.151 for the (M — TV, 7V)-soliton 
solutions obtained in Proposition |5.2[ Here we recall that 



C„ = / Vn+i{x,y,- ■ ■)dx, n=l,2,...,M, 



-oo 



where Vn+i = —dxSn{^d)lnT (see (3.14l). Note that for \y\ '3> the r-function is asymptotically 



equal to a single exponential function in each of the AI asymptotic sectors in the ccy-plane (see 



Figure 5.3). That is, in each asymptotic sector the r-function has the form 

T w ^{mi,...,mN)E{mi,...,mN) , 

and all other terms are exponentially small in comparison with this dominant exponential. Then we 
have 

N / oo \ 
Inr « E E Kn,tn]+C, 

where c is a constant determined by ^{mi, . . . ,mj^) and the /c-parameters. We now calculate C„ 
using Proposition |5.2| for y ^ and ?/ <C 0. 

For y 3> 0, there are iV -I- 1 asymptotic regions, and the boundary of adjacent regions is given 



by the line [e„,j„] for n = 1,....,N. Recall from Proposition 5.1 that the dominant exponential 
terms in adjacent regions differ by only one phase. We now integrate along a horizontal line 
which passes through those asymptotic regions, by using the above expression for In r and the fact 
that Sn{—d) — — ^9t„ + (h.o.d.). Since the contribution of the exponentially small terms vanish as 
y oo, we then obtain 

1 ^ 

Cn = -E {kl-kl) , n=l,2,..., 



=1 
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where each term in the sum is the contribution to the integral from the pair of dominant exponentials 
in adjacent asymptotic regions for y ^ 0. 

For y ^ 0, we follow the similar arguments as above and integrate w„+i along a line across the 
M — N+1 asymptotic regions separated by the lines [im , gm] ior m — 1, . . . , M — N . We then obtain 
another expression for the integral C„, namely 

1 M-N 

Cn = - E K-K) ' n = l,2,... . 

Since C„ do not depend on y, the two expressions obtained for y 3> and y <C must be the 
same, that is, we obtain the relations, 

1 N 1 M-N 

n ^' n s=i 

Rearranging the terms, we have for n — 1, 2, . . . , M 

M N M-N 

E fc" = E + E K 

j—1 r—1 s—1 

N M-N 

- E kl + E fc; 

r— 1 s — 1 



Then Lemma 5.1 implies that there is a one-to-one correspondence between two sets {fce^, fcg^} and 
{kj^^kij\. Furthermore, since {ke^,kgj\ is a set of distinct elements, then so is {kj^^,ki_^}. This 
proves that tt is a bijection. □ 



Note that the permutation tt defined by (5.6 1 has no fixed point because 7r(e„) = j„ > e„, n = 
1, . . . ,N and 'n{gm) — im < 9mi rn — 1, . . . ,M — N. Such permutations are called derangements. 
Moreover, tt have exactly N excedances defined as follows: an element I € [M] is an excedance of tt 



summarize the results of Proposition 5.2 and Corollary |5.1| as follows 



if it{1) > I. The excedance set of tt in (5.6) is the set of pivot indices {ei, 62, . . . , e^v}- We can now 



Theorem 5.4. Let A be an N x M irreducible matrix which gives a coordinate of a non-negative cell 



of the Grassmannian Gr(N,M). Then the tau-function (5.1) associated with the A-matrix generates 



an {M — N , N) -soliton solutions. The total M asymptotic line-solitons associated with each of these 



solutions induce a pairing map tt defined by {5.6). Moreover, tt is a derangement of the index set 



[AI] with N excedances given by the pivot indices {ei, 62, ... , ejv} of the A-matrix in RREF. 

The derangements tt S Sm are represented by linear chord diagrams with the arrows above the 
line pointing from e„ to j„ for n = 1^2, . . . ,N , while arrows below the line point from g,„ to i„i for 



TO = 1, 2, . . . , A/ — iV. Figure 5.4 illustrates the time evolution of a (3, 3)-soliton solution. The chord 



diagram shows all asymptotic line-solitons for y — > ±00. 

6. (2, 2)-SOLITON SOLUTIONS 

In this section we discuss all soliton solutions of the KP equation generated by the 2x4 irreducible 
A-matrices with nonnegative minors. Theorem |5 . 4| implies that each of the soliton solutions consists 
of two asymptotic line-solitons as y ^ ±00, that is, we have (2, 2)-soliton solutions. We outline 
below the classification scheme for the (2, 2)-soliton solutions. First we note that there are only two 
types of irreducible matrices for the size 2x4, i.e. 

1 -c -d\ /I a -c 

Q I a h ) "^^'^ 1^0 1 6 

The condition for nonnegative minors implies that the constants a, h, c and d must be non-negative. 
In the first case, one can easily see that ad = is impossible because then ^(3,4) < from the 
irreducibility. Then there are 5 cases for nonnegative minors with ad ^ 0, which are (1) ad — bc> 0, 
(2) ad — be — 0, (3) b — 0,c 0, (4) c = 0, 6 7^ 0, and (5) b — c = 0. For the second case, since 
ab ^ due to irreducibility, we have only two cases: (1) c — 0, (2) c 7^ 0. Thus we have total seven 
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Figure 5.4. The time evolution of a (3, 3)-soliton solution. The permutation of 
this solution is tt = (451263). The fc-parameters are chosen as (A;i, A;2, . . . , fee) = 
(— 1, — ^, 0, |, 1, |). The dominant exponential for a; <C is £'(1,2,5), and each 
dominant exponential is obtained through the derangement representing the so- 
lution, e.g. after crossing [5, 6]-soliton in the clockwise direction, the dominant 
exponential becomes £^(1,2,6). That is, [5,6]-sohton is given by the balance of 
those two exponentials. 



A = 



different types of ^-matrices, and each A-matrix gives a different (2, 2)-soliton solution. We now 
summarize the results for all seven cases with 2x4 irreducible and totally non-negative A-matrices: 

(a) TT = (3412): This case gives T-type 2-soliton solution, and the asymptotic line-solitons are 
[1,3]- and [2,4]-types for \y\ — > oo. This type of solution was first obtained as the solution 
of the Toda lattice hierarchy [4], and this is why we call it "T-type" (see also ^iTj). The 
A-matrix is given by 

'I -c -d> 
^0 1 a b 

where a,b,c,d > are free parameters with ad ^ be > 0. This is the generic solution on 
the maximum dimensional cell (four free parameters) of Gr(2,4), and contains all possible 
line-solitons in the process of interaction, i.e. [i, j]-solitons for any {i,j} pairs (see the next 
section for details). 

(b) TT — (4312): The asymptotic line-solitons are given by 

(i) [1,4]- and [2, 3]-solitons in y > 

(ii) [1,3]- and [2, 4]-solitons in y < 0. 
The ^-matrix is given by 

'1 ~b -c 



A 



1 a 



where a, fe, c > are free parameters. Note that two line-solitons for y <C arc the same types 
as in the T-type solitons. We also observe all possible line-solitons including intermediate 
ones, but not for all t unlike the T-type. 
(c) TT = (3421): The asymptotic line-solitons are given by 

(i) [1,3]- and [2, 4]-solitons in ?; > 

(ii) [1,4]- and [2, 3]-solitons in y < 0. 
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The >l-matrix is given by 

^ _ A -c 



^0 1 a b 

where a,b,c > axe positive free parameters. This solution can be considered as a dual of 
the previous case (b), that is, two sets of line-solitons for y ^ and y are exchanged. 



The example discussed in Example 5.2 corresponds to this solution (see Figure 5.2 1 
(d) TT — (2413): The asymptotic line-solitons are given by 

(i) [1, 2]- and [2, 4]-sohtons in y > 

(ii) [1,3]- and [3, 4]-solitons in y ^ 0. 
The A-matrix is given by 

'l -c -d 



A 



Ola b 



where a,b,c,d > with ^(3, 4) — ad — be = 0. 
(e) TT = (3142): The asymptotic line-solitons are given by 

(i) [1,3]- and [3, 4]-solitons in y ^ 

(ii) [1,2]- and [2, 4]-solitons in y < 0. 
The A-matrix is given by 

^ / 1 a — 



^0015 

where a, 6, c > 0. This solution is dual to the previous one (d) in the sense that the missing 
minors are switched by ^(3,4) ^ <^(1,2). This case will be further discussed in the next 
subsection to describe a connection with the Mach reflection in shallow water waves. 

(f) TT — (4321): This case gives P-type 2-soliton solution, and the asymptotic line-solitons are 
[1,4]- and [2,3]-types for \y\ oo. This type of solutions fits better with the physical 
assumption for the derivation of the KP equation, i.e. a quasi-two dimensionality with weak 
y-dependence. This is why we call it "P-type" (see [T7j). The situation is similar to the 
KdV case, for example, two solitons must have different amplitudes, A[i_4] > ^[2^3]. The 
A-matrix is given by 

'10 -fe^ 
^0 1 a 

(g) TT = (2143): This case gives 0-type 2-soliton solution, and the asymptotic line-solitons are 
[1,2]- and [3,4]-types for \y\ — > 00. The letter "O" for this type is due to the fact that 
this solution was originally found to describe the two-soliton solution, see for example |llj . 
(see also [TT]). Interaction properties for solitons with equal amplitude has been discussed 
using this solution. However this solution becomes singular, when those solitons are almost 
parallel to the y-axis; in this case, the order of the /c-parameters is reversed, that is, k2 > 
leading to £^(2,3) < 0). This is contrary to the assumption of the quasi-two dimensionality 
for the KP equation [52]. We will discuss this issue in the next section. The A-matrix is 
given by 

'1 a 0^ 



A 



^ '^0 15 

In Figure |6.1[ we list the chord diagrams for all those seven cases. One should note that any 
derangement of ^4 with exactly two excedances should be one of the graphs. This uniqueness in the 
general case has been used to count the number of totally non-negative Grassmann cells [501 US • 
We also mention that these seven types of soliton solutions have been found by applying recursive 
binary Darboux transformations (see Appendix in [5]). 

Now let us describe the details of some of the (2, 2)-soliton solutions, which will be important 
for an application of those solutions to shallow water problem discussed in the next Section. In 
particular, we explain how the A-matrix uniquely determines the structure of the corresponding 
soliton solution such as the location of the solitons and their phase shifts. 
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Figure 6.1. The chord diagrams for seven different types of (2, 2)-sohton solutions. 
Each diagram corresponds to a totally non-negative Grassmannian cell in Gr(2,4). 



6.1. O-type soliton solutions. This is the original two-soliton solution, and the solutions corre- 
spond to the chord diagram of tt = (2143). A solution of this type consists of two full line-solitons 
of [1,2] and [3,4] (see Figure 6.2 1. Note here that they have phase shifts due to their collision. Let 
us describe explicitly the structure of the solution of this type: The r-function defined in (5.1) for 
this case is given by 

T = E{1,3) + bE{lA) + aE{2, 3) -I- abE{2, 4) , 

where a,b > are the free parameters given in the A-matrix listed in the previous section. As we 
will show that those two parameters can be used to fix the locations of those solitons, that is, they 
are determined by the asymptotic data of the solution for large \y\. 

For the later application of the solution, we assume that [1, 2]-soliton has a "positive" y-component 
in the wave-vector (i.e. tan ^[12] < 0), and [3, 4]-soliton has a "negative" y-component, (i.e. 
tan 4] > 0, see Figure [?]). (Recall that any line-soliton has a "negative" x-component in the 
wave-vector.) Then for the region with large positive x, we have [1, 2]-soliton in y > and [3,4]- 
soliton in y < 0. Those solitons are obtained by the balance between two exponential terms in the 
T-function: 

For [1, 2]-soliton in a; > (and y ^ 0), we have the dominant balance between E{1, 4) and E{2, 4). 
Then the r-function can be written in the following form, 

T f=ibE{l,A) + abE{2,A) 

= 26e*''+5(«i+»2) ^Qgj^ 1 (^Q^ -02 + 6*^2 



12J ' 



which leads to the [1, 2]-soliton solution in the region near di w 62 for large x > 0, 
u = 2^2 Inr « ^{k2 - kifsech" ^ {0^ -62 + 0+) ■ 



Here the phase shift 9^2 (+ indicates a; > 0) is given by 



(6.1) 



3+ -I fc4-fcl , 

iTy = in z In a 

kA - fc. 



— ki 
kA — fco 



The parameter a determines the location of the [1, 2]-soliton. For example, one can choose a so that 
9^2 = Oi which implies that the [1, 2]-soliton is passing through the origin at t = 0. 

For [3, 4]-soliton in x > (and y <C 0), from the balance r w a£'(2,3) -I- abE{2,4), we have 



U « ^(fc4 - k3fsech^^i93 - 04 + ^3^4): 



SOLITON SOLUTIONS OF THE KP EQUATION 



35 




t = -10 t = t=10 



Figure 6.2. The time evolution of an 0-type soliton solution. Each E{i,j) in- 
dicates the dominant exponential in that region. The parameters are chosen as 
(fci, A;2, fca, fc4) = (— |, — J, J, |), so that both line-sohtons have the same ampli- 
tude, ^[1.2] ~ ^[3-4] ~ 2' ^^'^ directions of the wave-vectors are given by 
*[3,4l=-*[i,2]=tan-i|w68.2°. 



with the phase shift, 



3+ _ 1„ ^3 - ^2 , , . ,k3- k2 



(6.2) 61+4 = In - In b i.e. b = -^e'"'''^ . 

Thus the parameter b determines the location of [3, 4]-soliton, and one can choose appropriate b so 
that 9^^ — 0. Thus the parameters in the ^-matrix can be determined by the asymptotic data of 
the phase shifts 9^^ for x and \y\ » 0. 

Now calculating the phase shift for the [3, 4]-soliton in a; < 0, one can find the total phase shift 
along the line-soliton: From the balance between the terms with £'(1,3) and 6i?(l,4), and following 
the previous arguments, we obtain 

It is interesting to note that the sum of the phase shifts for those solitons are conserved in the 
y-direction for each A-matrix, 

(6.3) 9+2 + 0u = 0+, + 9^2 = In - ln(a6) . 

We measure the total phase shift for each soliton in the reference of the soliton in a; <C 0, that 
is, the total phase shift is defined by ^34 = 9^^ — 6*34 (recall that the phase shift is generated by 
interacting with other soliton, and the effect of the interaction propagates only in the "positive" 
x-direction) . We then obtain the well-known formula for the phase shift 6*34 (see for example '^^), 

{k4 - ki){k3 - k2) 

6^34 = In ■ 



{ks - ki){k4 - k2) 



Notice that this does not depend on the A-matrix. One can also see from (6.3 1 that the total phase 
shift of [1, 2]-soliton is the same as that of the [3, 4]-soliton, i.e. 

O12 = O12 ^ ^12 = ^34 = In Ao , 



where we note 



^ _ {kz - k2){ki - ki) _ {k2 - ki){ki - k^) ^ ^ 

(fc4 - fc2)(fc3 - fcl) (^4 - ^2)(fc3 - fcl) 



This implies that 9i2 = ^34 < 0, and each [i, j]-soliton shifts in x with 
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Figure 6.3. 0-type interaction for two equal amplitude solitons. The parameters 
fc,'s are taken by (fci, fcz, /cs, ^4) = (-1 - 10""*, -10"'', 10"'', 1 + 10"''), which give 
^0 = ^[1,2] = ^[3,4] = 5 and = tan [3,4] = -tan^-ji^a] = 1 + 2 x lO"* (i.e. 
^[3.4] ~ 45.0057). The constants a, b in the ^-matrix are chosen so that the center 
of interaction point is located at the origin (see (6.4|), and Mmax — u{0, 0, 0) « 1.96 
(about four times larger than each soliton amplitude Aq). 



The negative phase shifts Ax[i 2] < and Aa;[3 4] < indicate an attractive force in the interaction 

and if the amplitudes are the same, i.e. k2 — ki = — k^, then their phase shifts are the same 

1 



Aa 



[1:2] 



= ^2^[3,4] = 



fe,_fc3>934 < 0. Figure 
which have the same amplitude, ^[1.21 = A 



6.3 



[3,4] 



[1,2] 



^[3.4] 



illustrates an 0-type interaction of two solitons 
"Aq, and are symmetric with respect to the y-axis. 



0- 



0-type soliton solution has a steady X-shape with phase shifts in both line- solitons. One can also 
find the formula of the maximum amplitude which occurs at the center of intersection point (center 
of the X-shape): We place the soliton solution so that the origin (0, 0) is the center of the X-shape 
at i = 0. This implies that the sum of the phase shifts vanishes for each soliton, i.e. 



(6.4) 



9+ 
'12 



9+ ■ 
^34 



= ln(^i-^^#^-21na = 0, 



^34 



In 



{ki - fc2)(fc3 - ^2) 

(fca - fci)(fc3 - k2) 



21n6= 0. 



(^4 - fci)(fc4 - k2) 

Those determine a and b in the A-matrix. Then the r-function becomes 



(fc3 - fci) (^E^E; + E2Ei + ^K^{E^Ei + E2E^)) 
[E^+^/K^Ei 



El 

E- 



12 



^Ao £^3 + E, 



34 ■ 



EE cosh 6+ 



V £^34 

'Aq coshO 



) + -E'12 (^^£'34 + -^34 



where a = (3 implies that a and (3 give the same solution u, and Ef, and O* are defined by 



Ef, exp(±- 



e± _ [(0^ „ 02) ± (03 - 04)] . 
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It is then not difficult to see that u{x,y,0) attains the maximum at the origin, and we get the 
maximum value Wmax '■= u(0, 0, 0), which is 




k2){k3 - k4) 



(6.5) 



l + \/Ao 

(See also [Sl|32].) Since < Aq < 1, we have the bound 



^[1,2] + ^[3,4] < ""max < (^^^^[1.2] + \J ^[3,4] 



It is also interesting to note that the formula Aq has critical cases at the values fci = fc2 or 
^2 = fcs or k^ = k^. For the case with fci = fc2 or fca = fc4 (i.e. Aq = 1), one can see that one of the 
line-soliton becomes small, and the limit consists of just one-soliton solution. On the other hand, for 
the case k2 = fca (i.e. Aq = 0), the r-function has only three terms, which corresponds to a solution 
showing a Y-shape interaction (i.e. the phase shift becomes infinity and the middle portion of the 
interaction stretches to infinity). This limit has been discussed in as a resonant interaction of 

three waves to make Y-shape soliton. This limit gives a critical angle between those solitons which 
can be found as follows: First let us express each kj parameter in terms of the amplitude and the 
slope, 



fci,2 = ^ (tan«'[i,2] T y^2yl[i,2]) , 
^3,4 = ^ (tan [3,4] T y2yl[3,4] 



where the angle ^[ij] is measured in the counterclockwise direction from the y-axis (see Figure 6.3 1 
In particular, we have 



tan^'[i,2] -^2A[i,2] +2fc2, tan^'[3,4] = ^2A[3,4] + 2A:3. 

Without loss of generality, let us consider the special case when both solitons are of equal amplitude 
and symmetric with respect the y-axis i.e., ^[1.2] — ^[3,4] — ^0 and ^'[3,4] = — ^'[1,2] = ^'o > 0- This 
corresponds to setting fci — —k^ and ^2 = —k^. Then, for fixed amplitude Aq, the angle ^'o has a 
lower bound given by 

tan = ^2^0 + 2^3 > \/2Aq := tan *c ■ 
The lower bound is achieved in the limit k2 — k^ — 0, and the critical angle ^'c is given by 

(6.6) ^ c = taiT^ y/2AQ . 

Note that this function is monotone increasing in Aq, and for Aq = ^ the critical angle occurs 



when the two solitons intersect perpendicularly. In this symmetric case, from (6.5), the maximum 
amplitude given at the center of interaction can be calculated as 

4^0 _ . , 2Ao 



(6.7) U,nax = — 7^, with Ao = 1 - 2 ,T, ■ 

1 + v^o tan^ Wq 



Thus, at the critical angle — (i-e., Aq — 0), we have u,nax = 4Ao (see also 1^ 

We also note that there are three other soliton solutions of (2, 2)-type which have the same limit. 
These are the cases marked by tt = (3412), (2413) and (3142), and their chord diagrams degenerate 
to the same graph of figure eight in the limit. We will discuss these cases in the following sections 
Since (2413)- and (3142)-types are dual each other, we only discuss the (3142)-type. 

One should note that if we use the form of the 0-type solution even beyond the critical angle, i.e. 
^3 < ^2, then the solution becomes singular (note that the sign of E{2, 3) changes). In earlier works, 
this was considered to be an obstacle for using the KP equation to describe an interaction of two line- 
solitons with a smaller angle. On the contrary, the KP equation should give a better approximation 
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t = -20 t = t=10 



Figure 6.4. The time evolution of a (3142)-type soliton solution. The parameters 
are chosen as (fci, /c2, ^3, ^4) = ("§> ~5j 51 §)• The amplitudes of those solitons are 
^[1.3] = ^[2,4] — 2 and ^[1,2] — ^[3.4] — The directions of the wave- vectors are 
^[2,4] = ~^[i,3] = 45°. The critical angle is given by = tan"^ V^A) ~ 63.4°, 
which also gives ^'c = ^'[3,4] = —^[1.2]- 



to describe oblique interactions of solitons with smaller angles. Thus one should expect to have 
explicit solutions of the KP equation describing such phenomena. It turns out that the new types 
of (2, 2)-soliton solutions discussed above can indeed serve as good models for describing line-soliton 
interactions of solitons with small angles. We will show in the next section how these solutions are 
related to the Mach reflections discussed in [2^ . 

6.2. (3142)-type soliton solutions. We consider a solution of this type which consists of two line- 
solitons for large positive x and two other line-solitons for large negative x. We then assume that 
the slopes of two solitons in each region have opposite signs, i.e. one in ?/ > and other in y < 



(see Figure 6.4 1. The line-solitons for the (3142)-type solution can be determined from the balance 
between two appropriate exponential terms in its r-function which has the form, 

T = E{1, 3) + bE{l,A) + aE{2, 3) + abE{2, 4) + cE{3, 4) . 

The solution contains three free parameters a, b and c, which can be used to determine the locations 
of three (out of four) asymptotic line-solitons (e.g. two in a; 3> and one in a; <C 0). Thus, the 
parameters are completely determined from the asymptotic data on large \y\. 

Let us first consider the line-solitons in x ^ 0: There are two line-solitons which are [1, 3]-soliton 
in y ^ and [2, 4]-soliton in y ^ 0. The [1, 3]-soliton is obtained by the balance between the 
exponential terms &i?(l,4) and cE{3,A), and the [2, 4]-soliton is by the balance between aE{2,3) 
and ci?(3,4). Consequently, the phase shifts of [1,3]- and [2, 4] -solitons for x ^ are given by 



Choosing appropriate values for the ratios a/c and b/c, one can set 9^^ = 0^ = 0, thus fixing the 
locations of those line-solitons as passing through the origin at t = 0. Notice that there is one free 
parameter left after this. 

Now we consider the line-solitons in x <C 0: They are [3, 4]-soliton in y ^ and [1, 2]-soliton in 
J/ <C 0. The phase shifts are given respectively by 

6^., = In — — In b, OZ, = In — — In a 

The four phase shifts and the three free parameters a, 5, c are related by 

kl — kl n- _ krt — fcl n- fc-J — fci „+ „- 
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s<1 s=1 s>1 



Figure 6.5. A (3142)-type soliton solution with the s-parameter. The k- 
parameters are taken as (/ci, fc2, fcs, A;4) — (—2,-1,1,2). The parameters in the 
A-matrix are chosen as a = = |s and c — 3s, so that [1,3]- and [2, 4]-solitons 



meet at the origin, i.e. 6^^ — 6*^ ~ (see (6.10 1). Then at s = 1, all the solitons 



meet at the origin, i.e. the s-parameter shifts [1,2]- and [3, 4]-solitons. 



with the conservation of total phase shifts for y ±oo (as in the case of 0-type, see (6.3)) 



We then define the parameter representing this conservation, 

(6.9) s exp (-0+ - 6^,) , 
which leads to 

ks-ki g+ ks-ki g+ ks-ki 

(6.10) 0,= -, i—se^'^, b=- — sei-^, c=i ^ s. 

fes — k2 K4 — ki K4 — K3 

The s-parameter represents the relative locations of the intersection point of the [1,3]- and [3,4]- 
solitons with the x-axis, in particular, 9^^ + 9^^ = when s = 1. Thus the parameters a, b and c are 
related to the locations of [1, 3]-soliton (with 9^^), of [2, 4]-soliton (with ^?^), and the intersection 
point of [1,3]- and [3, 4] -solitons (with s). 

Now we consider an example, in which the [1,3]- and [2, 4]-solitons have the same amplitude 
(^[1,3] = ^[2,4] = ^0) they are symmetric with respect to the x-axis (^'[2.4] — ^^'[1,3] = ^o)- 
Then in terms of the /c-parameters, we have 

^3 — fcl = ^4 - 

Also the symmetry of the wave-vectors, i.e. ^'[2,4] 

^2 + ^4 = — (fci 
This implies that we have 

(6.11) ki = -ki > 0, fcg -fca > 0. 

The angle ^'o takes the value in (0, ^'c), where the critical angle is given by the condition k2 ~ k^~ 0, 
i.e. 

^ tan^i ^2^0 < |. 

Notice that this formula is the same as that of the 0-type soliton solution (see ( |6.6[ )). In this sense, 
these two cases have the same limit but from opposite directions namely, ^'o — > ^'c from above for 
the 0-type, while ^'o from below for the (3142)-type. This fact will be important for the 

initial value problem discussed in Section [7] 




From (6.11 1, one can easily deduce the following facts for [1, 2]- and [3, 4]-solitons in x < 0: 
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E(3,4) 



t--10 t = t=10 



Figure 6.6. The time evolution of a P-type soliton solution. The parameters are 
chosen as (fci, fc2, fcs, A:4) = (— |,— The amplitudes are ^[1,4] = ^,^[2.3] = 
|. The directions of the wave-vectors are given by ^'[2.3] = — ^'[1,4] = tan^^ ^ w 
26.57°. The parameters in the A-matrix are chosen as a = | and 6 = 3, so that the 
center of the interaction point is located at the origin at t = 0. 



(a) Those solitons have the same amplitude, i.e. 



1 



A[1^2] = ^[3,4] T^(^4 - h) 



Thus, if the [1, 3]- and [2, 4] -solitons in x > are close to the y-axis (i.e. a small ^'o), then 
the amplitudes of the solitons in a; < arc small; whereas at the critical angle ^0 = the 
solitons [1,2] and [3,4] in a: < take the maximum amplitude ^[1,2] — ^[3.4] — ^o- 
(b) The directions of the wave- vectors for the [1,2] and [3, 4] -solitons are also symmetric, i.e. 



tanvj/ 



[3,4] 



tan^fri 21 = k3 + ki. 



Moreover, the symmetry (6.11 1 implies that tan [3 4] = ^4 — ^2 = ^j2Ap^^ = ^/2Aq, so 



*[3,4] = = tan" 



/2Ao. 



Thus the directions of the wave- vectors for the [1, 2] and [3, 4]-solitons in a; < depend only 
on the amplitude of the solitons in x > but not on their directions (i.e., angle of their 
V-shape). 

Let us choose the parameters in the A-matrix for the (3142)-soliton solution appropriately, so 
that at i = all the solitons intersect at the origin (see Figure 6.4 1. Then for i > 0, the resonant 
interaction between [1,3]- and [3, 4]-solitons (as well as [2,4]- and [1, 2]-solitons) generates an inter- 
mediate line-soliton (called "stem" soliton) which is [1,4] soliton. The amplitude of this soliton is 
given by 

, 2 



(6.12) 



^1,4] = ^(fc4-fcl)' = 2fc2= ^ 



2Aq +tan^'o^ 

^I'ci the amplitude takes the maximum ^[1,4] = 4^o (see 



Note here that at the critical angle ^'g 

For t < 0, the resonant interaction between [1,3]- and [1, 2]-solitons (as well as [2,4]- and [3,4]- 
solitons) generates an intermediate line-soliton of [2, 3] -soliton. The amplitude of [2, 3]-soliton is 
given by 

1., . .0 1 ' ■ ^2 



A 



[2,3] 



'2Aq - tan 



Because of the symmetry (6.11 ), both [1,4]- and [2, 3]-solitons are parallel to the y-axis, i.e. tan [14] 
tan [2,3] = 0. 
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6.3. P-type soliton solutions. This type of solution consists of two line-solitons, [1,4] and [2,3]. 
We again assume that [1, 4]-soliton has positive slope and [2, 3]-soliton has negative slope (see Figure 



6.6). 



The phase shifts of [1,4]- and [2, 3]-soliton in a; ^ are given by 



'14 



In — - — In 6, 



, ^4 - ^2 , 

in ; in a. 



^4 — ^3 k4 — ^3 

Those parameters a and b in the ^-matrix can be chosen, so that at i = those solitons meet at 
the origin. 

The phase shifts of those solitons in a; ^ are 



^14 = In 7^ T' ~ ^' 

K4 - K2 

The total phase shift of [1, 4]-soliton is given by 



= in ; In a. 



714 



(fc2 - ki){ki - fca) 

Again this value is the same as the shift for [2, 3]-soliton, i.e. 623 = 6'^3 — 6*^3 = 014. Because of the 
order fci < A:2 < ^3 < A;4, the phase shift is always positive, i.e. 

1 



(6.13) 



Ax[i,4l = ^^^14 > 0. 



This situation is similar to the case of the KdV solitons. Notice in particular that the total phase 
shift has the opposite sign to that of the 0-type soliton solution ^ , which indicates a repulsive force 
in the interaction. 

As in the case of 0-type, one can find the amplitude at the center of interaction: We calculate 
the sums of Of^ and 6*^3, and set 

(fc2 - fci)(fc3 - fci) 
= Ti r\7T, - 2 In 6 = 0, 

(fc4 - fc2)(fc4 - K3) 

Q+ .a- 1 (^2 - ki){ki -k^) 

ij„ + c^23 = ill T, r^, - 2 In a = 0. 

^ ^ {ks - ki){k4 - ks) 

Then following the same procedure as in the case of 0-type, we get the amplitude at the center of 
interaction which is located at the point ofx^y = t = 0, 



'14 



'14 



m(0,0,0) 

(6.14) 

where Ap is given by 



I {{h - kif + {ks - k^y) - ] , "^ (^4 - h)ik3 - fcz) 



/At 



A 



[1,4] 



A 



[2,3] 



Ap = 



1 + VAp 

(fc2 - fcl)(fc4 - ^3) 

(fca - ki){k4 - k2) 



^[1,4] ^[2,3] , 



Ifwesetx= y 



then we have 



l~{x + y) 



k4 — ki 

l-VA^ ^ ^/{l-x){^-y)- ^ ^ 

l + ya^ " v/(l-x)(l-y) + " (v/(l-x)(l-j/) + V^)2 
Using the inequality (equivalent to ^/xy < ^{x + y)), 



V(l-a;)(l-2/) + V^< 1. 



we have 



/A^ 



>l-{x + y) 



A 



2,3] 



fc4 — fci 



^[1,4] 
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Figure 6.7. Two-soliton solution with T-type interaction. The parameters are 
chosen as (fci, fe, fca, fc4) = a.nd s = 1 (i.e. no phase shifts). Two 

other elements b, c in the A-matrix are chosen, so that two solitons intersect at the 
origin when r = 1. Also notice that T-type soliton consists of (2413)-soliton in 
a; < and (3142)-sohton in a; > 0. 



Then (6.141 gives the following bound of the interaction amplitude, 

2 



(\/^[l,4] - \/^[2,3]) < U(0,0,0) < A[i^4] - A[2,3] 



As a result, we find that the amplitude at the center of the amplitude is smaller than the larger 
soliton-amplitude ^[1,4] due to the P-type interaction which acts as a repulsive force (see also the 



phase shift (6.13|). 



6.4. T-type soliton solutions. There are four parameters in the A-matrix for T-type soliton 
solution consisting of two asymptotic line-solitons, [1,3] and [2,4]. Here we explain that those 
parameters give the information of the locations of those line-solitons, the phase shift and on-set of 



the opening of a box (see Figure 6.7). Thus three of those four parameters are determined by the 
asymptotic data on large \y\, and we need an internal data for the other one. 

Following the arguments in the previous section, one can find the phase shifts of the line-solitons 
of [1, 3] and [2, 4]: For [1, 3]-soliton in a; > (and y ^ 0), the phase shift is calculated as 

^13 = 1^ 1 ; 1^ IT ' 

where D — ad — be — ^(3,4). For the same soliton in a: < (and y <C 0), we have 
So the total phase shift is given by 

" " (fc2-fcl)(fc4-fc3) be 

Thus the shift depends on the A-matrix unlike the cases of O- and P- types, and it can take any 
value. 

For [2, 4]-soliton in a; > (and y -^0), we have 

.^, = ln^-ln^. 
ki~ k^ c 

and for the same one in a; < (and y ^ 0), we have 

^24 = In T T" ~ ^• 

fc4 — fci 
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Note that the total phase shift 624 = — O24 is the same as that for [1, 3]-soliton, i.e. the phase 
conservation along the y-axis ^13 + ^24 = ^r3 + ^24 holds. Then as in (6.9 1 for the case of (3142)-type, 
we define the s-parameter, 

(6.15) s:=exp(-0+ -^^24), 

which represents the intersection point of [1,3]- and [2, 4]-soliton. With the s-parameter, we have 

(6.16) ^=1 •se^^ c=- — se24, D^- — s. 

ki — ki kz — K2 ^4 — ^3 

Namely, the three parameters 6, c and D = ad — be determine the locations and the phase shift (i.e. 
the intersection point of [1,3]- and [2, 4]-solitons). One other parameter is then related to an on-set 
of a box at the intersection point (see Figure 6.7). 

In order to characterize this parameter, let us consider the intermediate solitons of [1,4] and [2, 3]. 
First note that for t ^ 0, [1, 4]-soliton appears as the dominant balance between E{1, 2) and E{2, 4). 
This implies that the r-function can be written by 

T « E{1, 2) + dE{2, 4) = cosh ^(6*1 - 6*4 + 61+ ) , 

where the super index means t > 0. Then the phase shift 9^^ is given by 

0+ =ln^^-lnd. 

K4 - K2 

Similarly one can get the phase shift 9^^ for t <C as 

ks-ki ^ D 



14 = ln- In—. 

kA — k-i a 



Now consider the sum of 9f^, 



n+,n-_ 1, (fc2-fcl)(fc3-fci) dD 
^" + ^"-^"(fc4^fc2)(fc4-fc3)" V- 

Also, for the [2, 3]-soliton, one can get 

(fc3 - fcl)(fc4 - ^3) d 

Now we introduce a parameter r in the form, 

a k^ — k2 



(6.17) 

so that we have 



d fcs — fci ' 



+9Y4 = ln-^, 9+ + 92, = - In (rs) . 



Suppose that at t = 0, [1, 3]- and [2, 4]-solitons in a; > are placed so that they meet at the origin, 
that is, we choose 9^, = 6*^ = 0. Also if there is no phase shifts for those solitons, i.e. s = 1. then 
the sums become 

9t, + 9^^ = lnr = -{9+, + 92,). 
This implies that at i = if r = 1, then the T-type soliton solution has an exact shape of "X" 
without any opening of a box at the intersection point on the origin. Moreover, at i = if r > 1, 
then [1, 4]-soliton appears in a; > and [2, 3]-soliton in a; < 0; whereas if < r < 1, then [1,4]- 
soliton appears in a; < and [2, 3] -soliton in a; > 0. Figure 6.7 illustrates those cases with s = 1. 
The parameter r determines the exponential term that is dominant in the region inside the box. 
When r < 1, E(2, 4) is the dominant exponential term, and when r > 1 the dominant exponential is 
E{1, 3). One should note that the parameter r cannot be determined by the asymptotic data, that 
is, r is considered as an "internal" parameter. 
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Figure 7.1. Initial data with V-shape wave. Each hne of the V-shape is locally a 
line-soliton solution, u{x, y, 0) = Aq sech^ ^^Ao{x + coy) with cq = tan(±^'o). We 
set those line-solitons to meet at the origin. In our numerical simulation, we reverse 
the time, i.e. t —t, so that the solitons are propagate in the positive x-direction, 
and small radiations propagate in the negative x-direction. 



7. Numerical simulation with V-shape initial waves 

In this section, we present some numerical simulations of the KP equation with obliquely incident 
initial waves related to a physical situation (see for examples [291 [33l US] ) • In these simulations, we 
change the sign of the time t, so that now the solitons move in the positive a;-direction (this seems 
to be a convention in the area of fluid mechanics). A complete numerical study of the initial value 
problem will be reported in a future communication |16j . 

We consider the initial data given in the shape of "V", 

(7.1) u{x,y,0) = Aq sech^y^i^o {x - co\y\) , 

where cq = tan'i'o > with the angle measured from the y-axis (see Figure 7.1 ). Note here that two 
semi-infinite line-solitons are propagating toward each other, so that they interact strongly at the 
corner of the V-shape. At the boundaries y = ±ymax of the numerical domain, those line-solitons 
are patched to the exact one-soliton solutions given by 



U(x, ij/inax, t) = Aq SCch^ ^ ^Aq {x =F Co?/max " l^t) 



with u = jCq + ^Aq. Notice that these solitons correspond to the exact one-soliton solution of the 
KdV equation. The numerical simulations are based on a spectral method with window-technique 
similar to the method used in [331 (the details will be reported in [16 ). The V-shape initial wave 
was first considered by M. Oikawa and H. Tsuji (see for example [29j i33]) in order to understand 
the generation of freak (or rogue) waves. They noticed generations of different types of asymptotic 
solutions, and found the resonant interactions which create localized high amplitude waves. In this 
paper, we present the results for the cases corresponding to = 2 and two different angles, ^'i 
and with < < ^'2. where the critical angle is given by = tan^^ V^A^ 63.4°. Then 
we explain these results in terms of certain (2, 2)-soliton solutions discussed in the previous section, 
and in particular, we describe the connection with the Mach reflection discussed in [25]. 

The main idea here is to consider the V-shape initial wave as the part of some (2, 2)-soliton 
solution listed in the previous section. In order to identify those soliton solutions from the V-shape 
initial wave form, let us first denote them as [ii, ji]-soliton for y ^ and [12, J2] -soliton for y <^ 0. 
Then using the relations, kj — k^ = \/2Aq and kj + fc,; = tan^E'oi for [i, j]-soliton, we have 



(7.2) 



%2 = o (*^^^ *o + 2) , = 5 (tan ^-0 - 2) , 



% = -;r(tan*o -2), = -i(tan«'o-l-2). 
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Because of the symmetry in the initial shape, we also have kj^ — —k^-^ and fc^^ = ^^ji- Moreover, at 
the critical angle, tan^fc — ^/2Aq — 2, we have ki^ — kj^ — 0. We also identify ki^ as the smallest 
parameter and kj^ as the largest one, so that depending on the angle ^'o, we obtain the following 
ordering in the fc-parameters: 

For < *o < *c (i-e. < tan^o < 2), we have 

ki-^ < ki.^ < < kj-^ < kj^ , 

implying that the corresponding chords of the [ii, ji]- and the [^2, j2]-solitons overlap. This means 
that the two solitons can be identified as part of either the (3412)-type (T-type) or the (3142)-type 



solution (see Figure 6.11. That is, [1,3] chord appears on the upper side of the diagram, and [2,4] 
chord on the lower side. As will be shown below, the numerical simulation suggests that the solution 
converges asymptotically to a (3142)-type soliton solution. This situation is related to the Mach 
reflection discussed in [HI [T^] . 

For < ^'o < f (i-e. 2 < tan^'o), we have 

In this case, the corresponding chords are separated implying that the two solitons form part of 
either (2413)- or (2143)-type (0-type) solution. That is, [1, 2]- and [3,4]-chords appear on the upper 
and lower sides of the chord diagram, respectively. The numerical simulation seems to indicate that 
the solution converges asymptotically to an 0-type solution, and this corresponds to the regular 
reflection [221IT2]. 

We can thus predict the following types of the asymptotic solutions depending on the values ^'q: 

(a) If the angle is ^'o = 0, then the initial wave is just one-soliton solution. This soliton is 
parallel to the y-axis, so fci = — ^2 showing no y-dependency. 

(b) If the angle is in the range < < = tan'^ \/2Aq, then the asymptotic solution is the 
(3142)-type sohton solution (not T-type) 

(c) If the angle satisfy ^'c < ^'o < f ; then the asymptotic solution is the (2143)-type soliton 
solution (i.e. 0-type, and not (2413)-type). 

In the next two subsections, we present the results for the cases (b) and (c) (the case (a) is 
obvious). 



7.1. Numerical simulation for '^^ > Vfc- In Figure [7T2| the upper three figures illustrate the 

numerical simulations of the initial value problem for the V-shape initial wave with Aq = 2 and 
= tan^^(3) w 71.57°. Here the critical angle is ^'c = tan~^(2) « 63.4°, and we observe 



asymptotically an 0-type soliton solution. The corresponding /c-parameters are obtained from (7.2 1, 
i.e. 

(^1, ^2; ^3) ^4) — (^1' 5; §)• 

Since two solitons [l,2]-and [3, 4]-solitons for a; > in the initial data are [1,2]- and [3, 4]-solitons, 
which meet at the origin, i.e. = ^34 = 0> take the A-matrix in the form, 

, /I I 



where we have used (6.1 1 and (6.2 1, i.e. a = T^E^ — | ^^^d b = — i. The corresponding 

0-type exact soliton solution is illustrated in the lower figures in Figure ]7?2l We also confirm that 
the numerical solution is indeed converging to this exact solution (details of these results will be 
reported in |16j). 

7.2. Numerical simulation for < ^c- In Figure |7.3| the upper figures illustrate the result 



of a simulation with the initial V-shape wave (7.11 with Aq ~ 2 and = 45°. The angle is 



now less than the critical angle ^'c. The asymptotic solution is expected to be of (3142)-type whose 
parameters are determined as follows: 
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Figure 7.2. Numerical simulation with the V-shape initial wave and an exact 
solution corresponding to the asymptotic solution. The parameters are given by 
Ao = 2 and ^-q = tan-i(3) = 71.6° > = tan-i(2) « 63.4°. The solution 
converges asymptotically to the 0-type solution shown in the three figures in the 
second line. The parameters are given by (fci, ^2, k^, ^4) = (— |, — |, |, |). 



(i) The fc-parameters can be obtained from (7.2 1 above 



(fcl,fc2, fcs, ^4) = (-§, -5, 5, §)■ 

The angles of the line-solitons of [1, 2]- and [3, 4]-types in x <C are then given by 
*[3,4] = -*[i,2] = = tan-i(2) » 63.4°. 

(ii) The elements of the A-matrix may be obtained by the condition that the intersection point 
at t = is located at the origin for the initial solitons identified as [1,3]- and [2,4]-types. 



Then from the phase shift formulas (6.8 1, we choose f = 3 and ^ = 1, so that [1,3]- and 
[2, 4]-solitons meet at the origin at t = 0. We also expect that all four solitons (two original 
in x > and two other solitons in a; < 0) may be considered to intersect at the origin at 



t — 0. This assumption leads to s = 1 in (|6.10|), and we obtain 

A- 



12 0-2 



1 I 



However this choice of the s-parameter may not be correct, since the generation of the line- 
solitons in X < affects the interaction timing at t (it seems that a correct choice is 
s < 1, and this will be discussed further in |.16j). 



The lower figures in Figure 7.3 illustrate the corresponding exact solution with this A-matrix. 
One can observe that the asymptotic solution with the V-shape initial wave with = 45° < ^'c 
converges to the corresponding (3142)-type soliton solution. An explicit error analysis for the results 
will be also reported in [IB] . 

7.3. The Mach reflection. In [22 , J. Miles considered an oblique interaction of two line-solitons 
using 0-type solutions. He observed that resonance occurs at the critical angle VPc, and when the 
angle is smaller than ^fc, the 0-type solution becomes singular (recall that at the critical angle 



SOLITON SOLUTIONS OF THE KP EQUATION 



47 




-20 20 -20 20 -20 20 




-20 20 -20 20 -20 20 



T=0 T=2.5 T=5 

Figure 7.3. Numerical simulation with the V-shape initial wave and an exact 
solution corresponding to the asymptotic solution. The parameters are given by 
Aq — 2 and ^'o — 45°. The intermediate wave has the amplitude about 4.5. This 
solution seems to converge asymptotically to the (3142)-type solution shown in the 
second line which is generated by the A-matrix with a = 2, 6 = | and c = 2 (i.e. 
all four line-solitons meet at the origin at f = 0). The parameters are given by 

(fci,fc2,fc3,fc4) = (-§ |, f). 



^'c, one of the exponential term in the r- function vanishes). He also noticed a similarity between 
this resonant interaction and the Mach reflection found in a shock wave interaction. This may be 
illustrated by the left figure of Figure |7.4| where an incidence wave shown by the vertical line is 
propagating to the right, and it hits a rigid wall with the angle — ^'o measured counterclockwise from 
the axis perpendicular to the wall (see also [l2]). If the angle of the incidence wave (equivalently 
the inclination angle of the wall) is large, the reflected wave behind the incidence wave has the same 
angle ^Pqi a regular reflection occurs. However, if the angle is small, then an intermediate wave 
called the Mach stem appears as shown in Figure |7.4[ The critical angle for the Mach reflection is 
given by the angle "^c- The Mach stem, the incident wave and the reflected wave interact resonantly, 
and those three waves form a resonant triplet. The right one in Figure 7.4 illustrates an equivalent 
system of the wave propagation in the left figure (one should ignore the effect of viscosity on the 
wall, i.e. no boundary layer). At the point O, the initial wave has V-shape with the angle ^'o, which 
forms the initial data for the simulation. Then the numerical simulation describes the refiection of 
line-soliton with an inclined wall, and these results explain well the appearance of the Mach refiection 
in terms of the exact soliton solution of (3142)-type. 

The maximum amplitude for this problem occurs at the wall, called the maximum run-up, and 



this can be obtained by the formula (6.7 1 for the 0-type solution, i.e. ^'o > ^'c, and the formula 
(6.121 for (3142)-type solution. Figure 7.5 illustrates the maximum amplitude for various angles ^'q. 
Notice that at the critical angle ^'c ~ 63.4°, the amplitude takes four times higher than that of the 
incident wave (this figure was first obtained in [H]). One can also find the length of the Mach stem 
(i.e. the intermediate soliton of [l,4]-type) from (3142)-soliton solution, that is, for the point {x,y) 
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Figure 7.4. The Mach reflection. The left figure illustrates a semi-infinite line- 
soliton propagating parallel to the wall. The right figure is an equivalent system 
to the left one when we ignore the viscous effect on the wall. The resulting wave 
pattern is of (3142)-soliton solution. 




Figure 7.5. The maximum amplitude obtained from the 0-type soliton for 
< ^'c and the (3142)-type soliton for > ^'c. The critical angle is 
= tan^^ ^2^0 ~ 63.4°. The highest amplitude is obtained at the critical angle, 
and it is 4 time larger than the initial soliton. 



of the resonant interaction of the triplet, we have 

1 o 

a; - (tan^fc + tan^fo) t, 



1 



(tan^fc - tan^'o) t, 



where tan 

by 



/2Ao and the point O is assigned as (0, 0). The angle <i> in Figure 7.4 is then given 

2(tan^'c - tanvPn) 
tan$- ^ ' 



{tan^c + tan^'o)^ ' 

This formula might be useful to determine the s-parameter for (3142)-type soliton solution. 
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